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Two  new  computer  codes,  DYNA2CYL  and  DYNASPHERE , have  been  developed 
for  use  in  system-ge neiated  Q4P  (SGJMP)  and  internal  EMF  (IIMP)  problems. 

The  code  DYNA2CYL  solves  Maxwell's  equations  dynamically  and  self- 
consistently  in  two  dimensions  (r  and  z)  in  the  region  between  two  concen- 
tric cylinders  of  finite  length.  The  ctalculations  are  for  end-on  irradia- 
tion of  the  cylinders,  which  is  simulated  by  specified  emission  of  electrons 
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20.  ABSTRACT  (Continued) 


from  either  or  both  the  outer  and  inner  cylinders  into  the  cylindrical 
cavity.  The  two  cylinders  can  be  isolated  from  ore  another  or  connected  by 
an  arbitral y load.  The  outputs  of  the  code  are  fields  and  currents  on  the 
inner  ^nd  outer  cylinders,  the  potential  difference  between  the  two  cylinders, 
and  the  current  through  the  load  between  the  cylinders. 

The  code  DYNASPHERf!  solves  Maxwell's  equations  dynamically  and  self- 
consistently  in  the  region  between  two  concentric  spheres  for  specified 
electron  emissions  from  the  inner  sphere  into  the  cavity.  The  outputs  oi 
this  code  are  similar  to  those  from  DYNAcCYL. 

Parameter  studies  are  performed,  varying  certain  of  the  essential  input 
parameters  3uch  as  characteristic  pulse  time,  fluence,  dimensions,  and 
emitted  electron  energy. 

The  codes  have  been  further  exercised  to  compare  quasi-static  solutions 
with  the  complete  dynamic  solutions.  Ranges  of  validity  of  the  quasi-static 
approximation  are  defined. 
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1.  INTRODUCTION 


This  document  inscribes  continuing  progress  in  the  area  of  computer 
code  dove lopments , modeling,  and  applications  for  problems  in  internal  l.MT 
(.11! MP ) and  system-generated  F.MP  (SGEMP).  This  work  is  an  outgrowth  and 
extension  of  the  effort  on  the  quasi -static  TFPIFM  codes  (Ref,  1 and 
Appendix  F)  . The  primary  accomplishments  are: 

1.  The  addition  of  new  geometries  of  particular  interest  to 
system- rel ated  assessments, 

2 Conversion  of  the  codes  to  a fully  dynamic  treatment, 

3,  Quant i f i cat  ion  of  errors  introduced  by  the  quasi  static 
approx  i mat  i ons , 

4,  A parameter  exercise  of  the  dynamic  codes  to  determine 
variations  of  currents  and  fields  as  a function  of  re’.  »- 
vanf  input  parameters  of  interest. 

The  PYNA2CV1.  code  was  developed  to  treat  coaxial  c\  Under..  ot  finite 
length,  one  located  within  the  other.  Flectron  emission  can  occur  from 
the  su  ice  of  either  ot  the  cylinders,  and  the  two  cylinders  may  he  con- 
nectc  y an  arbitrary  load  i njie  dance . This  g^ome^ry  is  useful  for  model- 
ing ollite  SGFMP  as  well  as  IFMP  in  cavities  containing  objects  such 
as  equipment  boxes. 

The  PYNASPHFRF  code  is  similar  to  PYNA2CYI.  and  was  developed  to  treat 
the  geometry  consisting  of  concentric  spheres  with  electron  emission  from 
the  inner  sphere.  The  two  spheres  may  le  either  isolated  or  connected  by 
a conductor.  The  solution  technique  consists  of  a finite-difference 
approach  to  the  complete  set  of  Maxwell's  equations.  Self-consistent 
particle  motion  is  considered  whereby  the  electron  trajectories  art* 
altered  hv  the  electric  field.  The  double-spherical  geometry  is  particu- 
larly useful  in  the  investigation  of  satellite  SGF.M1  analysis  and  simula- 
t ion  studies . 

PracsAf  pip  tu 


The  complete  s*t  of  c ode  geomet  ries  now  available  is  shown  in  Table 
1 -l.  Those  *odt?s  have  boon  dove  1 oped  od  refined  over  a period  of  sever,’’! 
years  and  have  been  applied  in  the  analysis  of  c number  of  j.yste«B  of 
interest.  Each  of  the  geometries  has  been  useful  for  different  aspects 
of  these  studies,  and  the  results  have  contributed  significantly  to  the 
understanding  of  lEMVSGEMP. 


Table  1-1 
AVAILABLE  COOES 


Code 

Geoart  »v 

Quasi -St at  tc 
or  Dynamic 

App 1 i cat  ions 

DIODE 

Parallel  plate 

Static 

i:ast-rur.ning  code  for  pill- 
box geometries 

SPARC 

Parallel  plate 

Quasi -st at i c 

As  above  only  quasi -st at ic 

TEDIEM-R6 

Infinite  cylinder, 
si  tie -on  illci, 

Quasi  - st at  ic 

Long  cylinders  such  as  mis- 
s i le  syst  ems 

TEDIEM-  R: 

Finite  cylinder, 
end-on  l 1 lum. 

Quas i -st at i c 

Small  cavities,  equipment 
boxes,  satellite  equipment 
hays 

TF  HI  EM- PC 

Rod  oi  cable  over 

ground  plane 

QuiiS  1ST  (at  1 C 

Perturbations  provided  by 
cables  ir.  space- charge 
region 

DYNACYL 

Single  finite  c>  » . 
end-on  ilium. 

Pynairi  c 

As  in  TEMEM-R: 

DYNA2CYL 

Concenr.ric  finite 
double  cylinders 

Dvnafai  c 

Most  versatile;  perturba- 
tiot  s nroduced  by  objects 
within  equipment  bays  or 
boxes , :•  at e 11  i to  SGEMP 

DYNxSFHFRfc 

Concentric  sphere 

Dyr.aai  c 

Satellite  SGEMP 

ABORC* 

Arbitrary  bod> 
o f x vo  1 lit  i on  c ode 

Dvnami  c 

Same  as  DYNA2CYL 

Under  development 

The  new  dynamic  capability  of  the  DYNACYi.  code  was  ued  to  detr  mire 
the  range  of  appl icabi lity.  of  the  quasi-static  approximation  by  com;  uring 
results  ',ith  the  existing  quasi-static  TEDIEM-R"  code.  Defining  this  range 
is  important  because  this  approx  mat  ion  is  use'-  extensively  in  svstem- 
related  calculations,  and  because  the  quasi-static  calculations  are  quite 
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straight  forward  and  can  b*  performed  quickly  on  i ne  back  of  the  proverbial 
wvtlopt  to  obtain  quick,  bttttr  than  order -of- magnitude  estiMte^ 

Final ‘y,  a number  of  parameter  studios  have  been  performed  to  demon- 
strate* the  code  capabilities  and  to  illustrate  how  fields  and  currents  are 
affected  by  important  input  parameters.  There  are  several  basic  input  var- 
iables of  interest,  including: 

Pulse  tiae  history  (e.g.,  rise  tier  and  pulse  width). 

Electron  energy  (or  velocity). 

Object  site. 

Photon  fluence  (or  magnitude  of  electron  emission  current). 

It  is  necessary  cwlv  to  solve  problems  for  a limited  range  of  these  vari- 
ables, Scaling  laws  (Appendix  A)  can  be  applied  to  extend  the  results  of 
the  calculations  to  other  dimensions  and  time  histories,  allowing  a few 
studies  to  describe  a wide  range  of  parameters. 

The  remainder  of  this  report  is  organ i:ed  as  to l laws. 

Section  2 — Increased  Coaqiut  at  ional  Capab.lity 

He  ct  i on  .t  - C caput  e r Code  Ve  r i f i c at  i on 

Section  4 — Applications  and  Parameter  Studies 

Section  S — Comparisons  of  Quasi-Static  and  Hulls  Dynamic 

So  I ut I on s f o r E l e c t romagn ct  u:  field  Cal  cm  I a* i on s 
in  a Cylindrical  Cavity 


■ c nee  d t * * ^ed<  i , n s'  iv  a -s  i n g » ■.  v 1 \ ph\  v i . .i  i ^ i t u.tt  i f oc  »!<' ,! 
the  haste  capabilities  : t lorsrr  computer  cod  s.  According  1'  , two  new 
cades  have  been  generated  for  use  in  calculating  electron  motion  and  field 
generat  ion.  These  add  to  our  previous  repertoire  (Refs.  1.;  and  Appendices 
l ,1.1  in  that  they  allow  the  solut  ion  to  the  tul  1 set  of  Maxwell  ' s equations 
for  geometries  where  there  are  two  separate  bodies.  both  of  the  new  codes 
are  for  rotational  ly  syma»et  nc  situations,  and  therefore,  consider  two 
spatial  dimensions,  ()ne  of  the  codes  is  written  in  cylindrical  geometry 
and  the  other  m spherical  gcomelrv.  Beth  codes  require  as  inputs  the 
emission  characteristics  of  the  electrons  leaving  the  \arious  surfaces. 

1 h«  c -.hIc'  ,lS!  h \.  i i I . written  in  C \ 1 t ml  r i c o 1 geonw  tr\  , .in  extension 

of  the  nVN’UYl  code  (Appendix  l')  to  include  an  inner  can.  The  numerics  of 

the  calculations!  procedure  are  very  similar  for  the  two.  However,  in  the 
new  version,  the  inner  can  is  represented  by  a roc  ion  if  high  conductivit'  , 
which  is  tantamount  to  setting  the  electric  field  to  rero  m a region  of 

the  inner  space.  A region  ot  specified  conduct i vi t v aiong  the  avis  con- 

nects the  inner  and  outer  cans  .in J allows  the  representation  of  a load 
between  the  two  conductors. 

Hie  mam  difference  in  the  two  cylindrical  codes  is  > n the  sophisti- 
cation ot  the  inputs  am«  outputs,  further  complicating  the  aituation  is 
the  t.’-.t  that  now  emission  can  take  place  from  either  the  inner  or  the 
outer  can  or  t rom  both.  Particle  tracking  is  also  more  difficult  since 
then*  is  now  an  inner  region  which  will  stop  any  incident  electrons.  How- 
ever, a wide  variety  of  problems  of  practical  interest  can  be  treated  hv 
using  the  increased  capability. 

Two  typical  problems  of  interest  are  shown  schemat i cal ly  in  figures 
2-1  and  2-2  ;>hotons,  being  a typical  cause  of  electron  emission,  are 
shown  ioi  illustrative  purposes  only,  since  the  emission  characteristics 
of  the  electrons  are  required  as  input. 
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RT-C^lft 

figure  1 , Schematic  representation  of  geomct  rv  for  cylindrical  code 
Sf'.fMP  calculations,  missile  compartment  with  equipment  box 

RESISTIVE  LOAD 


PHOTONS 

RT-09319 

figure  ^-2.  Schematic  representation  of  geometry  for  cylindrical  code 
SQ-'MP  calculat  ions;  configuration  for  StTVT  experiment 
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Figure  2-1  shows  the  sort  of  situation  one  might  expect  to  find  when 
a structure  is  inside  some  container.  The  structure  might  be  a satellite 
suspended  in  a test  chamber  for  photon  exposure,  or  it  might  be  a piece  of 
hardware  in  a conpartment . In  the  former  case,  the  lc>ad  represents  a cable 
monitoring  the  response  of  various  systems,  and  in  the  latter  case,  the 
load  represents  wiring. 

Figure  2-2  shows  a situation  where  the  interior  cylinder  has  been  dis- 
torted into  a disk  and  placed  at  one  of  the  outer  cylinder.  This  configu- 
ration is  of  particular  interest  for  an  experimental  situation  where  the 
end  of  a cavity  may  be  isolated  from  the  rest  of  the  system  by  a load 
resistance  to  allow  monitoring  of  the  collected  current.  The  influence 
of  voltage  build-up  on  the  collector  (in  that  it  can  change  the  current 
collected)  is  thus  automatically  accounted  for. 

The  code  DVNASPHERE,  written  in  spherical  geometry,  is  an  extension 
of  the  TSPHERE  code  (Appendix  E) . The  mathematical  description  of  the  code 
is  given  in  Appendix  B.  In  this  case,  the  region  between  two  perfectly 
conducting  concentric  spheres  is  treated.  At  present,  emission  is  con- 
sidered from  the  inner  sphere  only,  and  the  capability  of  including  a load 
between  che  spheres  does  not  exist.  Thus,  in  certain  aspects,  the  code  does 
not  have  the  diversity  of  the  UYNA2CYL  code.  However,  as  mentioned  in 
Appendix  B,  the  coordinate  system  for  the  field  calculations  is  included  in 
the  program  in  a very  general  manner:  the  basic  equations  are  written  in 

general  orthogonal  coordinates.  In  practical  terms,  this  means  that  vari- 
ations in  zone  size  are  relatively  easily  incorporated.  This  is  in  contrast 
to  the  situation  in  DYNACYL  and  DYNA2CYL,  where  the  constant  zone  size  is 
built  more  deeply  into  the  numerical  differencing  methods. 

This  ability  to  have  varying  zone  sizes  has  been  used  to  allow  treat- 
ment of  problems  with  high  current  densities.  These  cases  are  character- 
ized by  tiie  fact  that  there  -re  steep  gradients  near  the  emitting  surface 
requiring  fine  grid  spacing  (small  radial  zone  size)  there,  whereas  this 
requirement  does  not  hold  far  from  the  emitting  surface. 

The  two  computer  codes  have  broadened  our  computational  capability 
in  several  respects.  With  the  ability  to  treat  two  cylinders  connected 
by  a load,  we  can  describe  many  situations  of  interest  for  both  prediction 
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purposes  and  experimental  purposes.  With  the  ability  to  consider  variable 
zoning  in  the  radial  direction,  we  can  describe  cases  where  high  current 
densities  are  of  interest.  In  what  follows,  results  of  various  applica- 
tions of  these  capabilities  are  shown. 


3.  COMPUTER  CODE  VERIFICATION 


Before  the  results  of  a computer  code  can  be  accepted,  and  prior  to 
application  of  a code  to  systems  studies,  the  code  must  be  tested.  These 
tests  usually  amount  to  comparison  of  results  with  analytical  solutions 
where  available,  and  conparison  with  other  code  results  for  regions  where 
both  codes  are  applicable. 

The  DYNA2CYL  code  is  a relatively  direct  extension  of  DYNACYL  as  far 
as  the  physics  involved  is  concerned.  Verification  of  the  new  code  is 
accomplished  by  comparing  the  results  with  the  old  code  where  possible. 
This  comparison  is  presented  below.  Results  are  also  presented  and  inter- 
preted showing  the  logical  consistency  of  the  code.  Finally,  DYNASPHERE 
results  are  conpared  with  those  of  two  other  codes.  One  of  these  is 
TSPHERE , which  differs  from  the  new  code  in  that  it  uses  a quasi -static 
solution  to  Maxwell's  equations.  The  second  of  these,  LFLUX,  is  valid 
only  for  low  fluences  but  solves  the  full  set  of  Maxwell's  equations. 

3,  1 DYNA2CYL  CODE  CHFCKO’JT 

The  DYNA2CYL  code  is  a modification  of  DYNACYL  in  which  an  inner  con- 
ducting cylinder  has  been  placed  within  the  simple  cylindrical  cavity 
treated  by  DYNACYL.  Numerous  checkouts  have  been  performed  to  ensure 
proper  code  operation.  Comparisons  with  DYNACYL  (empty  cavity)  have  been 
made  for  the  case  of  increasingly  smaller  inner  cylinders  to  ensure  that 
the  solutions  approach  the  proper  limit.  Comparisons  between  the  surface 
currents  on  the  inner  cylinder  have  been  made  with  currents  on  the  surface 
of  a spherical-shaped  object.  The  magnitudes  and  time  histories  (taking 
into  account  retardation)  agreed  well  for  low-  and  mediura-f luence  pulses 
having  a 20-nsec  FWHM.  Comparisons  were  not  made  for  high  fluence,  but 
good  agreement  is  not  expected  for  such  cases  since  dynamic  and  geometric 
effects  will  be  important  due  to  the  shortening  of  the  effective  pulse 
length  by  the  space  charge  limiting. 
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3.2  DYNA2CYL  DOUBLE- CYLINDER  IEMP  STUDY 


3.2.1  Introduct ion 

In  modeling  compartments  of  a missile  or  satellite  for  IEMP  calcula- 
tions, complex  geometries  are  reduced  to  simple  ones  such  as  spheres  or 
cylinders.  Often  a cavity  which  is  cylindrical  in  shape  may  contain  equip- 
ment boxes,  greatly  conpl icating  the  geometry.  These  equipment  boxes  may 
be  connnected  to  the  outer  cavity  walls  by  circuits  which  have  effective 
impedances  varying  from  essentially  short-circuit  values  to  open-circuit 
values.  Thus,  the  effect  of  these  objects  in  the  cavity  can  greatly  com- 
plicate the  IEMP  response  of  the  system.  With  the  introduction  of  the 
double -cylinder  IEMP  capability  of  the  DYNA2CYL  code  comes  the  opportunity 
to  make  predictions  of  the  effects  of  objects  in  a cavity  connected  by 
arbitrary  circuits  to  the  cavity  walls.  The  code  is  currently  limited  to 
ar.  interior  cylinder  anywhere  on  the  axis  of  the  outer  cylinder  and  con- 
nected by  a resistive  load  to  the  outer  cylinder  (s=*e  Section  2).  This 
geometry  and  circuitry  can  model  many  cases  of  interest  with  much  greater 
accuracy  than  was  formerly  available. 

Results  of  calculations  performed  with  DYNA2CYL,  preserved  in  the 
following  subsections,  wi li  show  the  effects  on  IEMP  of  placing  an  object 
such  as  an  equipment  box  inside  the  cavity. 

Results  are  presented  for  end-on  irradiation  for  three  cases: 

1.  An  empty  cylinder, 

2.  A cylinder  within  a cylinder  (open-circuit) 

3.  A cylinder  within  a cylinder  (short-circuit). 

Results  for  both  moderate  and  high  space-charge  limiting  are  presented. 
Moderate  space-charge  limiting  occurs  when  the  ratio  of  peak  transmitted 
current  reaching  the  rear  face  to  the  emitted  current  is  approximately  half 
the  value  ac  low-fluence  emission.  High  space-charge  limiting  is  the  des- 
ignation used  when  the  ratio  of  transmitted  to  emitted  current  is  less  than 
10%  of  the  value  at  low  levels. 

The  values  chosen  for  cavity  size,  emitted  electron  energy,  and  time 
history  are  typical  for  IEMP.  Although  emission  from  all  faces  occurs  in 
practical  problems,  forward  emission  from  the  outer  cylinder  face  only  is 
considered  here.  This  facilitates  the  interpretation  of  the  results  and 


is  not  unrealistic  when  *everse  emission  from  objects  within  the  cavity  is 
less  than  the  forward-emitted  currents  from  the  outer  container  walls. 

3.2.2  Problem  Definition 

The  basic  test  problem  is  a cylindrical  cavity  with  a length  and 

diameter  equal  to  40  cm.  Monoenergetic  electrons  with  a velocity  of 
8 

1 x 10  m/sec  (20  keV)  are  emitted  within  the  cavity,  with  a triangular 
pulse  chape.  A pulse  rise  time  is  chosen  to  be  approximately  10  times  the 
transit  time  for  light  across  the  cavity,  a situation  which  will  produce 
a quasi-static  response  at  moderately  low  fluences  and  dynamic  response 
at  highei  fluences.  Hence,  the  pulse  rise  time  is  chosen  to  be  approxi- 

2 

mately  10  nsec.  Peak  emission  current  density  levels  are  0.1  and  1 amp/cm  , 
corresponding  to  moderate  and  high  space-charge  limiting.  For  problems 
in  which  an  inne  cylinder  is  considered,  a cylinder  of  20  cm  diameter 
and  length  has  been  chosen. 

Fields  and  currents  are  predicted  throughout  the  cavity.  However,  the 
electric  field  at  the  emitting  face  and  the  magnetic  field  at  the  rear  sur- 
frce  are  of  particular  interest.  The  entire  problem  geometry  and  field 
locations  are  shown  in  Figure  3-1. 

3.3  RESULTS  OF  CALCULATIONS 
3.3.1  Moderate-Fluence  Results 

Results  for  an  open-  and  short-circuited  cylinder  within  a cylinder 
are  compared  to  results  obtained  with  an  empty  cylinder.  Of  primary  inter- 
est is  the  electric  field  at  the  emitting  surface  and  the  magnetic  field 
at  the  rear  surface  wall. 

There  exists  a printout  option  in  the  DYNA  codes  which  allows  time- 
dependent  tracking  of  electron  trajectories.  Electron  trajectories  for 
the  three  cases  identified  above  are  shown  in  Figure  3-2,  Electron  paths 
are  shown  for  times  near  the  peak  of  the  emitted  electron  pulse.  Because 
of  the  slow  variation  of  the  emitted  electron  pulse,  the  trajectories  shown 
corx^spond  to  quasi-steady-state. 

Brief  comments  about  the  cases  follow. 
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ELECTRON  ENERGY  = 20  keV 
PULSE  FWHM  = 10  nsec 
EMISSION  CURRENT  LEVELS  = 0.1,  1 amp/cm2 


Configuration  used  to  determine  the  effects  of  the 
inner  cylinder  with  the  DYNA2CYL  code 


Figure  3-2.  Particle  trajectories  for  cases  of  an  empty  cylinder,  open- 
circuit  inner  cylinder,  tnd  short-circuited  inner  cylinder 
for  low  current  density 
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a. 


a.  Empty  Cylinder 

A modicum  of  space-charge  limiting  tends  to  bend  electron  trajectories 
towards  the  side  of  the  can  and  reduce  current  transmitted  to  the  rear  face. 
In  a strictly  one-dimensional  calculation,  a fraction  of  the  predicted 
charge  would  be  returned  (incorrectly)  to  the  emitting  face. 

b.  Open-Circuit  Double  Cylinder 

Charge  striking  the  inner  cylinder  cannot  return  to  the  emitting  walls, 
and  therefore,  charge  continually  builds  up.  However,  sufficient  charge 
does  not  accumulate  on  the  inner  cylinder  in  the  time  frame  of  the  pulse 
(moderate  fluence  only)  to  severely  alter  the  electron  trajectories. 

c.  Short-Circuit  Double  Cylinder 

Short-circuiting  the  inner  cylinder  effectively  reduces  the  potential 
at  the  center  of  the  cavity,  thus  pulling  electrons  toward  the  center 
cylinder. 

The  electric  and  magnetic  fields  for  each  of  the  three  cases  are 
shown  in  Figures  3-3  and  3-4.  Comments  regarding  differences  in  the  fields 
Follow . 

Electric  Fields 

The  electric  fields  for  the  empty  cylinder  and  open-circuit  cylinder 
are  essentially  the  same  for  the  first  part  of  the  pulse.  The  same  amount 
of  charge  exists  in  the  cavity  regardless  of  whether  it  is  in  the  free 
space  or  contained  on  the  inner  conducting  cylinder.  At  late  times,  the 
charge  in  both  cases  tends  to  fall  as  electrons  leave  the  cavity.  However, 
the  field  remains  at  a significant  level  and  in  fact  will  reach  a steady- 
state  (non-zero)  value  at  extremely  late  times  for  the  open-circuit  cyl- 
inder case.  This  is,  of  course,  because  charge  is  trapped  on  the  isolated 
cylinder. 

For  the  short-circuit  cylinder,  the  electric  field  at  the  face 
increases  for  a time  approximately  equal  to  the  electron  flight  time  from 
the  emitting  face  to  the  inner  cylinder.  After  the  charge  strikes  the 
inner  cylinder,  it  quickly  leaves  the  cavity  (transit  time  for  light  is 
short  compared  to  electron  transit  time  in  this  problem),  and  the  electric 
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field  does  not  substantially  increase  thereafter,  although  there  is  some 
increase  due  to  charge  *illing  the  remainder  of  the  cavity. 

Magnetic  Fields 

The  magnetic  field  is  considerably  reduced  for  the  open-circuit  object 
within  the  cylinder  when  compared  to  the  empty  cylinder.  The  fields  are 
reduced  because  the  net  current  transmitted  to  the  rear  face  is  substan- 
tially reduced  by  the  shadowing  effect  cf  the  inner  cylinder.  On  the  other 
hand,  tne  magnetic  field  is  increased  by  almost  a factor  of  two  when  the 
center  cylinder  is  short-circuited  to  the  rear  face.  The  presence  of  the 
inner  cylinder  considerably  reduces  the  space-charge  limiting  and  causes 
more  current  to  flow. 

3.3.2  High-Fluence  Results 

Th  ■ same  problem  is  now  treated,  except  that  the  emission  currents 
aie  increased  by  an  order  of  magnitude  to  produce  a high  degree  of  space- 
charge  limiting.  Electron  trajectories  corresponding  to  the  open  and 
short-circuited  inner  cylinder  cases  are  shown  in  Figure  3-5. 

a * Open- Circuit  Cy linden 

Figure  3-S  shows  electron  trajectories  for  the  oper.-circui t case  for 
electrons  emitted  at  5 nsec,  or  half-way  into  the  pulse.  Charge  initially 
strides  the  inner  c>  Under  and  produces  extremely  high  fields  and  poten- 
tials, causing  the  electrons  tc  be  deflected  to  either  the  emitting  face 
or  the  side  wai  ’s . 

b . Short-Circuit  Cylinder 

The  gro’Uiding  of  the  inner  cylinder  permits  the  charge  striking  the 
cylinder  to  flow  to  ground.  The  corresponding  low  fields  produce  smaller 
perturbations  tr»  electron  trajectories. 

The  electric  and  Magnetic  fields  corresponding  to  these  two  cases  are 
shown  in  Figures  3-6  and  3-7,  and  are  discussed  briefly  below. 

F. lectric  Fields 

The  electric  fields  for  the  open-  and  short-circuit  cases  appear  to 
be  surprisingly  similar.  One  would  at  first  think  that  the  electric  field 
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Figure  3-5.  Particle  trajectories  for  cases  of  a shorted  inner 
cylinder  and  open-circuit  inner  cylinder  for  high 
current  density 
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Electric  field  at  center  of  emitting  face  of  exterior 
cylinder  for  two  cases,  ’.igh  current  density 
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Figure  3-7.  Magnetic  field  at  outer  wall  of  exterior  cylinder 
at  the  end  away  frori  the  emitting  surface  for  two 
cases,  high  current  density 
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for  th«  short-circuit  cylinder  would  So  sipifiemt  ly  tool  tor  than  that 
for  tho  open-circuit  cylinder.  ISe  fact  Is  that,  for  highly  United  con* 
ditions,  the  electric  field  at  the  emitting  surface  is  quite  insensitive 
to  the  conditions  in  the  remainder  of  the  or  ty.  The  field  distributions 
in  space  and  potential  differences  between  the  enitting  surface  and  points 
in  the  cavity  do,  however,  depend  strongly  on  the  details  in  the  cavity. 
Thus,  the  electric  field  at  the  enitting  surface  is  not  a good  indicator 
of  the  degree  of  space-charge  Uniting. 

The  effective  decrease  i.i  rise  tine  due  to  the  varly  onset  of  space- 
charge  Uniting  excites  resonant  cavity  modes  and  p reduces  an  oscillatory 
response  in  both  electric  and  magnetic  fields. 

Magnetic  Fields 

Because  of  the  high  degree  of  space-charge  Uniting  for  the  open- 
circuit  cylinder,  relatively  snail  amounts  of  current  are  transmitted 
across  the  cavity.  The  magnetic  fields  are  lower  than  for  the  short- 
circuit  cylinder. 

3.3.3  Summary  of  DYNA2CYL  Results 

In  summary.  the  fields  and  currents  within  a cavity  equipment  box, 
satellite,  etc.,  can  be  significantly  affected  by  the  presence  of  objects 
and  the  grounding  used  for  these  objects.  The  response  in  terms  of  mag- 
netic fields  and  transmitted  currents  can  he  significantly  increased  by 
the  reduction  of  the  space-charge  barrier.  On  the  other  hand,  highly 
space-charge- 1 inited  electric  fields  in  the  vicinity  of  emitting  surfaces 
(e.g,  , where  many  systems  cables  nay  h?  located)  are  insensitive  to  objects 
within  the  cavity.  There,  a fairly  accurate  description  of  fields  within 
the  cavity  may  be  possible  even  when  the  cavity  is  quite  irregular. 


3.4  DYNASPHERE  CODE  CHECK-OlTT 

In  this  section,  comparison  is  nade  between  DYNASPtERE  and  the  quasi- 
static code  TSPHERE  to  check  the  accuracy  of  self-consistent  routines,  in 
the  following  subsections,  comparison  is  made  with  the  semi -analytic  code 
LFLUX  {Ref.  3)  *t  low  flue nee  where  space-charge  limiting  is  not  a 
consideration. 
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The  configuration  used  for  comparison  is  shown  in  Figure  1-1,  in  which 
•loctrons  era  ooittod  from  tho  inner  half-sphere,  Surface  currants  and 
electric  fields  ora  computed  and  coaparad  in  check-out  problems,  Compari- 
sons ara  prasantad  between  tha  total  surface  currant  on  the  sphere  passing 
across  a plana  through  tha  sphere  at  90*  from  the  incident  photon  flux. 


Figure  3-S.  Configuration  for  comparisons  of  the  dynamic  DYNASPHERE 
code  with  the  quasi -static  TSPHERE  code  and  the  In- 
fluence LFLUX  code 


3.  «.  I Comparison  of  DYNASPHERE  and  TSPHERF.  for  SCI  Conditions 

Surface  currents  obtained  on  the  inner  sphere  of  a concentric  sphere 
geometry,  using  the  newly  developed  DYNASPHERE  code,  have  been  compared 
with  results  obtained  from  the  earlier  TSPHERE  code  under  both  non-space- 
charge-  limited  (NSCL)  and  space-charge- limited  (SCI)  conditions.  TSPHERE 
is  a concentric- sphere  SCElf*  code  which  is  the  forerunner  to  DYNASPHERE. 

A cosplete  description  of  TSPHERE  can  he  found  in  Reference  2.  TSPHERE 
and  DYNASPHERE  treat  similar  geometries,  hut  TSPHERE  employs  the  quasi - 
static  approximation  to  obtain  electric  fields  and  DYNASPHERE  alloys  the 
full  Maxwell's  equation  s^t  to  obtain  the  fields. 
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Results  obtained  with  the  two  codes  are  expected  to  agree  for  long- 
pulse-  length  low-fluence  problems  where  the  quasi-static  approximation  has 
been  shown  to  be  valid.  Comparisons  were  performed  for  a concentric  spher- 
ical region  with  an  inner  radius  of  2 meters  and  an  outer  radius  of  5 meters. 
A slowly  rising  pulse  with  a long  pulse  width  was  chosen  for  purposes  of 
comparing  with  the  quasi-static  code.  The  emission  current  pulse  is  trian- 
gular with  100-nsec  rise  and  fall  times.  The  pulse  rise  time  is  approxi- 
mately equal  to  the  electron  flight  time  across  the  cavity.  A spectrum  of 
electrons  with  an  average  energy  of  5 keV  was  emitted  outward  from  the  inner 
sphere  with  a level  proportional  to  cos  0,  falling  off  to  zero  at  the  side 

(0  = 90°)  of  the  inner  sphere.  Peak  emission  current  densities  of  0.42 
2 

and  4.2  amp/m  were  chosen,  corresponding  to  little  space-charge  limiting 
and  moderate  space-charge  limiting. 

Electrons  were  emitted  from  six  points  on  the  forward  half  of  the 
inner  sphere.  TSPHERE  employed  17  radial  and  six  angular  zones.  The  radial 
zones  were  finer  near  the  inner  sphere  and  coarser  near  the  outer  sphere. 

This  zoning  was  used  to  describe  the  steep  gradients  which  occur  near  the 
inner  sphere  under  high-fluence  conditions.  DYNASPHERE  employed  17  x 6 
zoning  grids  with  equal  spacing  everywhere  for  the  low-  and  medium-fluence 
cases.  Thirty  radial  and  six  angular  zones  equally  spaced  were  employed 
for  the  high-fluence  case.  The  emission  electron  energy  spectrum  was  broken 
into  12  bins.  Results  were  obtained  for  the  first  350  nsec  the  pulse. 

a.  Low-Fluence  Results 

Results  for  the  surface  currents  on  the  inner  sphere  at  90°  to  the 
incident  photon  direction  are  shown  in  Figure  3-9.  Agreement  between  the 
codes  for  the  low-fluence  case  seems  to  be  good.  This  agreement  is  expected 
because  of  the  long  pulse  length  and  lack  of  fields  affecting  particle  tra- 
jectories. Note  the  time  labelled  ttR/2c  on  the  graphs.  This  is  the  time 
required  for  light  to  travel  from  the  0=0  position  to  the  0 = 90°  posi- 
tion on  the  inner  sphere.  The  DYNASPHERE  surface  current  displays  a time 
lag  evidenced  by  its  values  being  somewhat  lower  than  TSPHERE  results  dur- 
ing this  time  interval.  TSPHERE  does  not  take  this  retardation  into  account, 
so  its  surface  currents  rise  more  steeply  at  the  beginning  of  the  pulse. 
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Figure  3-9.  Surface  currents  at  90°  versus  ti  ; for  DYNASPHERE  and  TSPHERE 
codes  applied  to  long-pulse-lengti  problem,  low  fluence 


Note  that  the  peak  of  the  surface  current  occurs  at  a time  later  than 
the  peak  of  the  emitted  current  pulse.  The  reason  for  this  is  that  the 
electrons  rust  move  out  almost  one  spherical  radius  before  they  induce  the 
maximum  current  on  the  sphere.  Thus,  the  lag  between  the  emission  current 
peak  and  the  surface  current  peak  is  approximately  equal  to  the  tine  required 
for  the  electron  to  move  out  to  a large  fraction  of  the  sphere's  radius. 

b.  Moderate- Fluence  Case 

At  higher  fluences  there  is  space-charge  limiting.  The  emitted  pulse 
of  current  begins  rising  until  a point  is  reached  at  which  the  currents 
pulled  back  to  the  sphere  are  approximately  equal  to  the  emitted  currents. 

The  currents  induced  on  the  sphere  do  not  substantially  increase  after 
this  point,  and  therefore,  the  effective  rise  time  of  the  driving  function 
(emitted  electrons)  is  decreased.  When  this  happens,  cavity  oscillations 
may  be  introduced.  This  is  the  case  here,  and  the  oscillations  are  shown 
in  Figure  3-10.  Note  that  both  the  quasi-static  TSPHERE  and  dynamic  DYNA- 
SPHERE show  the  decreased  pulse  rise  time  (compare  rise  times  in  Figures 
3-9  and  3-10),  but  that  the  DYNASPHERE  code  displays  an  oscillating  solu- 
tion characteristic  of  the  dynamic  response  of  the  cavity,  whereas  TSPHERE 
does  not.  The  theoretical  characteristic  period  of  oscillation  for  a double 
spherical  cavity  (radii  2 and  5 meters)  is  approximately  50  nsec,  a value 
which  compares  favorably  with  the  results  of  Figure  3-10. 

3.4.2  Comparison  of  DYNASPHERE  and  LFLUX  for 
Short-Pulse-Length  Problem,  NSCL 

The  comparison  of  DYNASPHERE  with  an  existing  SGEMP  code  reported  in 
the  previous  section  was  valuable  for  testing  the  new  code  for  consistency 
with  earlier  methods  for  long-pul^e- length  problems  and  space-charge- limited 
currents.  That  comparison  does  not  provide  a test  of  the  dynamic  capabil- 
ities of  the  new  code,  however.  To  obtain  this  check-out,  DYNASPHERE  was 
compared  with  the  existing  LFLUX  code  on  a moderately  short-pulse-length 
problem.  LFLUX  is  a concentric-sphere  SGEMP  code  in  whicn  electron  motion 
is  not  consistent  with  the  electric  and  magnetic  fields.  Therefore,  its 
results  are  valid  only  for  low-fluence  conditions. 


SURFACE  CURKEN 


Figure  3-10.  Surface  currents  at  90°  versus  time  for  DYNASPHERE  and  TSPHERE 
codes  applied  to  long-pulse- length  problem,  medium  fluence 


The  conparison  of  DYNASPHERE  with  LFLUX  was  made  on  a low-fluence  prob- 
lem in  which  the  pulse  width  was  compirable  to  the  L/c  time  of  the  cavity. 
The  L/c  tins  is  the  characteristic  time  for  light  to  traverse  the  cavity. 
Studies  made  with  the  IRT  DYNACYL  code  for  the  solution  of  Maxwell's  equa- 
tions in  finite  cylinders  have  shown  that  pulse  widths  comparable  to  the 
photon  flight  time  can  cause  large  effects  directly  attributable  to  the 
dynamic  terms  in  Maxwell's  equations  (see  Section  5).  This  test  problem 
should,  therefore,  be  a good  check-out  of  the  dynamic  capabilities  of  the 
DYNASPHERE  code. 

The  configuration  of  the  test  problem  is  illustrated  in  Figure  3-8  and 
in  the  inset  of  Figure  3-11.  The  inner  and  outer  sphere  radii  were  2.3  and 
15.2  meters,  respectively.  Electrons  of  7-keV  energy  were  emitted  straight 
out  from  the  forward  half  of  the  inner  sphere  with  a uniform  spatial  distri- 
bution. The  pulse  had  a 45-nsec  rise  time,  10-nsec  flat  top,  and  45-nsec 

fall  time.  The  photon  and  electron  flight  times  were  roughly  40  and  200 

2 

nsec,  respectively.  The  peak  emission  current  density  was  0.11  amp/m  , 
which  resulted  in  no  effect  of  the  fields  on  electron  trajectories. 

In  the  DYNASPHERE  runs,  the  region  between  the  spheres  was  broken  into 
45  radial  and  10  angular  zones  of  equal  spacing.  Electrons  were  emitted 
from  the  centers  of  the  angular  zones,  which  were  in  the  forward  direction. 
The  time  step  employed  was  5 nsec,  which  allowed  for  a reasonable  descrip- 
tion of  the  emission-electron  pulse  shape  without  causing  large  computer 
costs.  LFLUX  employed  approximately  120  radial  zones.  The  angular  depen- 
dence in  LFLUX  is  done  semi- ana lytic'.lly,  employing  Legendre  polynomials. 

In  this  case,  the  first  five  were  used,  which  was  equivalent  to  six  angular 
zones  in  the  direct-differencing  method.  The  tine  step  was  1.33  nsec.  The 
velocity  of  the  electrons  is  5 x 10  m/sec.  A quick  calculation  will  show 
that  both  the  DYNASPHERE  and  LFLUX  time  steps  were  consistent  with  their 
radial  zone  sizes  in  that  the  electrons  moved  approximately  one  zone  each 
time  step.  DYNASPHERE  zoning  is  somewhat  rougher  than  LFLUX  zoning  because 
its  capability  to  handle  the  entire  fluence  range  requires  particle- 
following to  obtain  currents  as  opposed  to  analytically  specifying  them. 

This  increased  capability  results  in  higher  computer  costs. 
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The  results  of  the  coef>«ri.*on  are  shown  in  Figure  3-11.  Surface  cur- 
rents at  90*  on  the  inner  sphere  obtained  from  the  two  codes  are  plotted 
versus  tine.  Agreement  is  seen  to  be  excellent.  The  minor  disagreements 
may  be  attributable  to  the  rougher  grid  sizes  employed  by  DYNASPHERE.  Both 
codes  show  roughly  the  same  values  for  the  peak  surface  current  and  also 
the  peak  values  in  the  oscillations  due  to  cavity  resonances.  Also,  the 
period  of  the  resonances  is  the  same  in  both  codes.  This  period  corresponds 
roughly  to  the  time  for  light  to  traverse  the  cavity,  which  is  the  period 
expected  for  the  oscillations. 
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Figure  3-11.  Comparison  of  surface  currents  at  90°  versus  time  from  the 
1ST  DYNASPHERE  code  and  the  Mission  Research  LFLUX  code. 
Medium-pulse- length  case  with  pulse  length  <<  electron 
flight  ti«  , low  fluence. 


4.  APPLICATIONS  AND  PARAMETER  STUDIES 


Dynamic,  t ime-dependent  codes  now  exist  to  treat  several  geometries 
of  interest  for  systems  and  phenomenology  studies.  These  codes  have  been 
exercised  over  a range  of  parameters  of  interest  to  investigate  the  char- 
acteristics of  the  solutions  and  sensitivities  to  the  various  input  param- 
eters. A number  of  results  from  these  parameter  studies  are  presented  here. 

The  primary  input  parameters  of  interest  are: 

1.  Emitted  electron  energy  (or  velocity)  spectrum, 

2.  Dimensions, 

3.  Pulse  time  history, 

4.  Emitted  current  levels  (or  photon  fluence) . 

Many  variations  are  possible,  including  spatial  distribution  of  emit- 
ted electrons,  geometry  variations,  emissions  from  various  surfaces,  time- 
dependent  emitted-electron  energy  spectra,  effects  of  time-phasing  electron 
emission  (i.e.,  angle  of  incidence  effects),  etc.  The  study  reported  here 
is  restricted  to  the  four  basic  input  variables  identified  above  because 
they  appear  to  be  the  most  fundamental  as  well  as  the  most  interesting. 

Also,  certain  of  the  parameter  variations  mentioned  above  (such  as  effects 
of  emitted  angular  distribution)  were  reported  previously  (Ref.  3). 

4.1  APPLICATIONS 

Two  fundamentally  different  types  of  problems  can  be  treated  by 
either  the  DYNACYL  or  DYNASPHERE  code.  One  consists  of  the  "outside 
problem,"  in  which  electrons  are  emitted  from  the  surface  of  the  inner 
object  towards  the  outside  surface.  The  inner  object  is  representative 
of  an  isolated  body  such  as  a satellite,  while  the  outer  object  can  be 
representative  of  a satellite  test  chamber.  The  external  response  of 
the  object  in  free  space  may  be  investigated  by  moving  the  outer  boundary 
out  to  a very  large  distance,  thereby  approximating  a surface  at  infinity. 
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The  second  type  of  problem  is  the  "inside  problem,"  consisting  of 
enission  from  the  outer  well  towards  the  inside  in  which  there  night  or 
night  not  be  an  inner  object.  If  present,  the  inner  object  night  or  night 
not  be  connected  to  the  outer  bouidary  wall.  The  outer  cavity  wall  is  rep- 
resentative of  an  equipnent  box,  an  electronics  box,  or  a uissile  cavity, 
while  the  inner  object  night  be  an  equipnent  box  or  other  structural 
element. 

The  outside  problen  described  above  has  become  known  as  the  systea- 
generated  electromagnetic  pulse  (SGEMP)  problem,  while  the  inside  problen 
has  retained  the  name  of  internal  electromagnetic  pulse  (IEW)  problen. 

The  nonenclature  is  admittedly  confusing,  and  the  distinction  between  the 
two  problems  diminishes  as  the  cavity  becomes  more  and  more  open  to  the 
outside. 

In  this  section,  we  consider  only  the  outside  problen  and  treat  the 
inside  problem  in  the  next  section.  The  DYNASPHERE  code  is  exercised 
' ing  a wide  range  of  variables. 

In  the  following  studies,  calculations  were  made  for  a 2-meter- radius 
* ere  with  electron  emission  occurring  from  one-half  the  sphere,  as  shown 
in  Figure  3-3.  The  electron  energy  distribution  is  approximately  Maxwel- 
lian with  an  average  velocity  of  0.1  to  0.2c  unless  otherwise  noted. 

Pul st  time  histories  are  triangular,  with  rise  times  equal  to  the  pulse 
wiv  s.  Electron  fluence  levels  are  variable  and  have  beer  adjusted  to 
correspond  to  either  low,  medium,  or  high  space-charge  limiting.  "low- 
-.luence"  solutions  correspond  to  emission  levels  where  there  is  no  space- 
char^  limiting  and  the  solution  is  essentially  linear.  All  charge  emit- 
ted from  the  inner  sphere  is  transmitted  to  the  outer  sphere.  "Medium 
fluence"  corresponds  to  the  situation  where  v.  modicum  of  limiting  occurs 
and  approximately  10%  of  the  emitted  charge  is  able  to  escape  from  the 
inner  sphere  to  the  outer  sphere.  "High  fluence"  corresponds  to  a great 
deal  of  limiting  whereby  only  1%  or  less  of  the  emitted  charge  escapes. 

In  all  cases,  the  response  parameter  of  interest  has  been  chosen  to 
be  the  surface  current  crossing  a plane  which  passes  through  the  sphere 
at  an  angle  0 as  shown  in  Figure  3-8. 
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Before  presenting  the  results  of  the  parameter  study,  a few  remark* 
about  scaling  laws  will  be  made  to  indicate  how  the  results  of  the  present 
study  can  be  generalized. 

4.2  SCALING  LAWS* 

A detailed  discussion  of  scaling  laws  is  given  in  Appendix  A.  A sum- 
mary of  these  results  follows. 

Solutions  can  be  obtained  directly  for  a whole  class  of  problems  from 
a single  solution  obtained  for  one  set  of  parameters.  Thus,  the  solution 
for  a single  set  of  parameters  has  a wider  range  of  applicability  than  for 
the  particulai  set  of  parameters  chosen.  Solutions  are  identical  if  the 
input  parameters  are  changed  by  the  following  factors. 

t"  = t/o  (time) 

RJ  = R/a  (dimension) 


I*'  = af 


v = v 


(f luence) 

(electron  velocity) 


The  solution  scales  if  the  time  and  dimensions  are  compressed  by  a, 
the  incident  f luence  (or  total  number  of  electrons)  is  increased  by  the 
same  factor  a,  and  the  electron  velocity  (energy)  is  unchanged.  When  the 
input  parameters  are  varied  in  this  first  fashion,  the  electromagnetic 
field  quantities  will  change  in  the  following  way. 

E'  * aE  (electric  field) 

H'  ■ aH  (magnetic  field  strength) 

V'  = V (potential) 

I'  = I (current) 

2 

J'  = a^J  (current  density) 

-> 

p'  = a^p  (charge  density) 


As  an  example,  consider  the  solution  to  a problem  obtained  with  a 
pulse  of  characteristic  time  x,  dimension  L,  peak  current  density  .1^,  and 
electron  velocity  v.  The  solutions  ob  lined  for  this  problem  are  denoted 
by  E,  H,  V,  I,  J,  and  p.  For  problems  in  which  the  pulse  characteristic 
time  is  shortened  by  a,  the  dimensions  made  smaller  by  a,  and  the  peak 

*The  work  described  in  Section  4.2  and  in  Appendix  A was  carried  out 
at  IRT,  supported  by  IRSD  funds. 


emission  current  density  increased  by  (one  factor  of  o for  the  increase 
in  fiuence  and  one  factor  for  the  decrease  in  pulse  width) , the  solutions 
will  be  oE,  oH,  !,  o^J,  and  a‘p. 

The  solutions  presented  in  the  following  sections  were  obtained  for  a 
specific  set  of  parameters  characteristic  of  those  found  in  certain  problems 
of  interest.  However,  these  results  can  be  directly  scaled  to  obtain  solu- 
tions for  many  other  ranges  of  parameters  of  interest. 

1 

4.3  PARAMETER  STUDY  OF  RISE  TIME 
4.3,1  General  Background 

A parameter  study  of  pulse  rise  time  has  been  performed  in  which  the 
rise  time  and  pulse  width  (e.g.,  primary  frequency  components)  are  varied 
through  a range  where  dynamic  or  resonant  effects  may  he  important.  Reso- 
nances are  usually  excited  when  the  excitation  pulse  has  significant  energy 
content  in  the  frequency  range  of  the  structural  resonances.  Generally, 
this  occurs  when  the  pulse  has  a rise  time  and  pulse  width  shorter  than  the 
transit  time  of  light  across  tne  object  of  dimension  L, 

tr  > L/C  . 

Higher- order  modes,  of  course,  can  he  excited  with  pulses  of  shorter  rise 
times  (higher- frequency  content'). 

The  lowest-order  mode  for  a sphere  of  radius  R is  approximately 


which  corresponds  to  a frequency  of  0,8c/2wR.  The  lowest-order  mode  for 
the  sphere  can  be  strongly  excited  with  a pulse  having  a characteristic 
time  in  the  range*  of 

0.34  0.27R 

Tr  a ~ f c * 

Thus,  a pulse  of  approximately  16-nsec  rise  tine  and  width  may  couple 
efficiently  to  the  2-meter-radius  sphere  considered  in  this  parameter 

■ i.^V. 

*\  . ! 'e  of  r^se  time  tf  *•  .s  a significant  amount  of  energy  close 
to  the  fi  .ui  rev  f * 0.34/tr. 
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The  previous  discussion  pen  tins  to  an  Isolated  sphere  in  free  space. 
However,  DYNASPHEWE  solutions  are  obtained  with  an  outer  spherical  boundary. 
This  boundary  will  not  perturb  the  solution  for  the  isolated  sphere  for  the 
period  of  tine  required  for  light  to  traverse  to  the  outer  wall  and  back  to 
the  inner  surface.  This  period  of  tine  is  designated  "clear  tiae"  and  is 
simply 

2(Rq  - R) 
tc c • 

where  R^  is  the  outer  boundary  wall  radius. 

Resonant  frequencies  for  concentric  spheres  are  shown  in  Figure  4-1 
as  a function  of  outer  spherical  radius  for  an  inner  sphere  radius  of  2 
aeters.  The  priaary  frequency  content  in  SO- * JO- , and  lS-nsec  pulses  are 
marked  on  the  vertical  scale  for  co«$>arison . 

Surface  currents  for  a 2-aeter  sphere  (outer  boundary  radius  of  10 
aeters)  as  a function  of  time  are  shown  in  Figure  4-2  for  various  angular 
slices  around  the  sphere.  The  pulse  width  and  rise  tiae  of  the  eaitted 
electron  pulse  is  SO  nsec.  It  is  evident  from  Figure  4-1  that  such  a pulse 
has  a frequency  content  below  the  lowest -order  aode  associated  with  the 
double-spherical  cavity  and,  therefore,  should  not  strongly  excit.1  the 
cavity.  Indeed,  the  solutions  shown  in  Figure  4-2  do  net  display  any 
oscillatory  behavior. 

We  now  point  out  several  characteristics  of  these  solutions.  Note 
that  th.  peak  in  surface  current  occurs  at  approximately  20  nsec  after  the 
peak  of  the  emitted  electron  pulse.  This  delay  time  cannot  be  attributed 
to  the  finite  speed  of  light  between  9-0  and  Q » 90",  which  is  only  on 
the  order  of  t^  (w/2)( R/c)  = 10  nsec.  Furthermore,  if  it  were  due  to 
light  delay,  the  surface  currents  would  peak  at  different  times  for  the 
different  angles. 

The  delay  is  caused  by  the  finite  tiae  required  for  the  electrons  to 
move  away  from  the  sphere  a sufficient  distance  to  induce  maximum  surface 
currents  on  the  sphere.  If  we  consider  a single  electron  moving  away  from 
a sphere,  we  find  thet  the  maximum  surface  current  is  induced  immediately 
after  the  electron  leaves  the  surface.  However,  if  we  consider  a stream 
of  electrons  moving  away  from  the  surface,  we  find  that  each  electron 
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Resonant  frequency  of  cavities  associated  with  the  TM  Modes  of 
concentric  spheres;  inner  radius  fixed  at  2 Meters 


Figure  4-2.  Surface  currents  at  various  locations  on  a 2 -Meter- radius  sphere 


contributes  significantly  to  the  surface  current  until  it  reaches  1/2  to 
S/4  of  a radius.  The  surface  current  continues  to  rise  as  electrons  are 
emitted  wtil  the  electron  atrea » reaches  1/2  to  S/4  of  a radius,  at  which 
point  the  surface  current  reaches  its  peak.  Thus,  the  delay  tine  between 
the  mission  peak  and  the  surface  current  peak  is  simply  the  tine  required 
for  the  electron  to  reach  a distance  of  approximately  1/2R. 

For  the  exalte  problem,  we  are  considering  this  delay  time  as 

• rife  • 20  ns'c  • 

Thus,  a delay  time  of  20  nsec  added  to  the  50-nscc  pulse  rise  time  results 
in  a peak  surface  current  at  approximately  70  nsec,  as  shown. 

A second  characteristic  of  interest  is  the  ratio  of  peak  emitted  cur- 
rent to  peak  surface  current,  he  define  the  peak  emitted  current  1 _ to 
be  the  sum  of  all  peak  surface  current  densities;  i.c., 

'«• " j-’p  rr,V  dJ  • 

where  t^  is  the  time  at  which  the  local  emission  current  density  is  peak. 
The  peak  surface  current  is  the  total  peak  current  flowing  across  a plane 
at  the  angle  where  the  maximum  current  flows;  i.e., 

I - (6  ) - IK  r (&  ,t  ) dl  . 
surf  nr  J surf  v m’  p' 

The  angle  6^  at  which  the  peak  surface  current  flows  is  dependent  tg>on  the 
assumed  emission  pattern.  For  half-sphere  emission  assumed  in  these  stud- 
ies, the  maximum  occurs  close  to  90*.  It  is  evident  in  Figure  4-2  that 
the  surface  currents  at  0 ■ 45*  and  155*  are  lower  than  the  current  at 
6-90*. 

The  ratio  of  peak  surface  current  to  peak  emission  current  is  depen- 
dent on  several  factors  including  total  emitted  charge,  pulse  rise  time. 

and  electron  velocity.  For  highly  space- charge- limited  problems,  the 
-4 

ratio  can  be  10  or  less.  Even  for  low-emission- level  situations  where 
no  charge  returns,  the  ratio  can  be  significantly  less  than  unity  for 
fast  pulses  and  slow  electrons.  The  ratio  can  approach  0.5  for  half- 
sphere  emission  when  the  pulse  rise  time  is  long  cohered  to  the  time 
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required  for  the  electrons  to  travel  to  1/2  to  3/4  of  the  spherical 
radius*  (see  Appendix  C) . 


For  the  present  example  in  which  there  is  no  space-charge  limiting,  the 

rise  time  would  have  to  be  greater  than  2 meters  divided  by  vg  = 0.1c,  or 

66  nsec,  for  the  ratio  I J I to  approach  1/2.  The  rise  time  is  50  nsec, 

sui'r  em  rr 

and  therefore,  the  ratio  should  be  somewhat  less  than  0.45.  The  actual 
ratio  is  ^3.85/14.5,  or  0.27. 


4.3.2  Rise  Time  Variation 

Surface  currents  at  0 » 90°  ar-.  shown  together  in  Figure  4-3  for 
pulse  rise  times  of  50,  30,  and  IS  nsec.  The  pulse  width  is  fixed  at 
50  nsec,  and  the  peak  emission  current  density  is  identical  for  all  three. 
We  note  that  the  rise  time  of  the  surface  current  decreases  with  decreas- 
ing photon  pulse  rise  times.  The  delay  time  between  the  peak  in  surface 
current  and  the  peak  in  emitted  current  remains  on  the  order  of  20  to  25 
nsec  for  all  three  cases. 

Rtsonance  begins  to  appear  as  the  pulse  rise  time  is  decreased  from 
50  to  15  nsec.  This  resonance  is  associated  with  the  presence  of  the 
outer  sphere  rather  than  an  isolated  sphere.  Note  the  slight  increase 
in  current  at  a time  slightly  later  than  the  t - required  for  the  wave 
to  travel  to  the  outer  wall  and  back  to  the  inner  sphere. 

Strong  resonances  are  not  excited  for  any  of  the  pulses  due  to  the 
slow  velocity  of  the  electrons.  It  must  be  remembered  that  the  peak  sur- 
face current  does  not  occur  before  the  electrons  move  out  approximately 
1/2  to  3/4  of  a radius  regardless  of  pulse  rise  time.  For  slow  electrons 
(v  = 0.1c),  this  minimum  time  is  greater  than  the  time  required  to  excite 
resonances.  Thus,  we  see  evidence  that  the  quasi-static  solution  pre- 
vails when  the  electron  transit  time  across  a characteristic  dimension  of 
interest  (^R)  is  slow  compared  to  the  transit  time  for  light,  even  when  the 
pulse  rise  time  is  fast.  This  result  is  demonstrated  in  the  next  section. 
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Figure  4-3.  Surface  current  at  90°  location  on  a 2-meter- 

radius  sphere;  constant  pulse  width  and  fluence 
various  rise  times 


4.4  ELECTRON  VELOCITY  VARIATIONS 


In  the  previous  section,  it  was  stated  that  the  electron  velocity  is 
a controlling  factor  in  determining  resonance  effects.  In  this  section, 
we  explicitly  demonstrate  this  fact  by  comparing  two  sample  calculations 
for  which  the  mean  electron  speed  has  values  of  ^0.15c  and  ^0.9c,  with 
all  other  parameters  the  same. 

A 20-nsec  pulse  width  and  rise  time  were  chosen.  A low  emission  cur- 
rent density  was  chosen  so  that  neither  problem  would  be  space-charge- 
limited.  The  only  difference  in  the  input  parameters  to  the  two  problems 
is  the  velocity  distribution  of  the  photo-emitted  electrons. 

An  outer  boundary  radius  of  22  meters  was  chosen  so  that  enough  clear 
time  (^120  nsec)  is  available  to  observe  resonances  of  an  isolated  sphere. 
The  results  of  this  exercise  are  shown  in  Figure  4-4.  It  is  evident  that 
the  surface  current  corresponding  to  the  lower-velocity  electrons  does  not 
exhibit  an  oscillatory  response,  whereas  the  current  corresponding  to  the 
fast  electrons  is  clearly  oscillatory.  Note  also  that  the  delay  between 
the  peak  of  the  surface  current  and  the  peak  of  the  pulse  is  diminished 
for  the  high-velocity  electron  problem;  the  electrons  reach  i/2  to  3/4  of 
a radius  sooner.  The  surface  current  responses  in  Figure  >~4  are  normal- 
ised to  the  same  peak  value  to  emphasize  the  difference  in  the  oscillatory 
behavior.  The  actual  peak  value  for  the  "fast"  electrons  exceeds  that  of 
the  "slow"  electrons  by  about  a factor  of  three. 

The  preceding  results  would  be  much  more  pronounced  for  an  emitter 
with  a geometry  associated  with  high-Q  resonances.  It  is  well  known  that 
sphere  resonances  are  highly  radiation-damped,  whereas,  for  example, 
diunbbe  11  -shaped  objects  resonate  with  a much  lower  frequency,  in  propor- 
tion to  their  characteristic  dimensions,  and  ring  down  relatively  slowly 
compared  to  the  sphere  case.  It  is  not  our  purpose  here  to  demonstrate 
extreme  resonance  effects,  but  rather  to  present  the  fundamental  physics 
of  SGEMP- induced  surface  current  resonances. 

An  important  point  worth  repeating  is  that  the  characteristic  time 
for  determining  the  effective  spectral  content  of  an  SGEMP  excitation  is 
related  to  the  time  required  for  the  photo-emitted  electrons  to  move  out 
an  effective  distance  from  the  system,  in  addition  to  the  rise  time  of 
the  photon  pulse. 


Surface  current  (normalized)  at  90°  for  fast  and  slow  photo-emitted 
electrons;  fast  and  slow  curves  are  normalized  to  the  same  peak  value. 


4.S  FLUENCE  VARIATIONS 


All  the  previous  parameter  variations  have  been  performed  fbr  low- 
fluence  electron  emission  where  virtually  no  space-charge  Uniting  occurs. 
In  the  next  series  of  problems,  we  begin  with  a low  fluence  and  progres- 
sively increase  the  fluence  while  holding  all  other  input  parameters  con- 
stant. A pulse  width  and  rise  tine  of  20  nser-  have  been  chosen.  The 
inner  radius  is  2 neters  and  the  outer  radius  10  meters,  corresponding 

to  a clear  tine  of  approximately  50  nsec.  Eaission  currents  vary  froa 
2 4 2 

a low  of  3.2  aap/a  up  to  a high  of  3.2  x 10  aap/a  . The  corresponding 
surface  currents  are  shown  in  Figure  4-5.  Oscillatory  behavior  after  50 
nsec  is  associated  with  electromagnetic  waves  between  the  inner  and  outer 
walls,  while  oscillatory  behavior  before  50  nsec,  exhibited  at  the  higher 
fiuences,  is  associated  with  the  inner  sphere  alone. 

We  note  several  interesting  aspects  of  these  solutions  in  the  follow- 
ing subsections. 

4.5.1  Peak  Surface  Currents 

The  peak  surface  current  as  a function  of  emission  current  is  shown 

in  Figure  4-6.  The  onset  of  space-charge  limiting  occurs  at  an  emission 

2 

current  of  several  anp/m  , and  a further  increase  of  4 orders  of  magnitude 
in  emission  current  results  in  an  increase  of  a factor  of  50  in  surface 
current.  Thus,  space-charge  limiting  is  very  effective  in  limiting  the 
replacement  currents  flowing  on  an  object. 

4.5.2  Rise  Times  of  Surface  Currents 

It  is  interesting  to  note  the  effective  decrease  in  response  rise 
time  as  the  fluence  is  increased.  The  time  at  which  the  peak  surface 
current  occurs  is  shown  as  a function  of  fluence  in  Figure  4-7.  Note 
that  the  effective  rise  time  of  the  response  can  be  r.ignificantly  shorter 
than  the  rise  time  of  the  emitted  electron  pulse.  This  shortening  of  the 
surface  current  rise  time  occurs  because  the  space-charge- limited  elec- 
trons do  not  travel  far  from  the  surface  before  they  are  turned  around 
and  returned  to  the  surface.  The  response  is,  therefore,  significantly 
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Figure  4-6,  Peak  surface  current  on  a 2 -meter- radius  sphere 


Figure  4-7.  Time  of  occurrence  of  peak  surface  current 
on  a 2-meter- radius  sphere 
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smaller  Chan  it  would  he  in  the  absence  of  space-charge  Halting,  and  it 
is  also  significantly  faster  (with  aore  high-frequency  spectral  content). 

4.6  DIMENSIONAL  VARIATIONS 

All  previous  solutions  have  been  obtained  with  an  inner  sphere  of  2 
aeters.  Solutions  for  other  dimensions  can  be  obtained  directly  by  using 
the  scaling  laws  described  in  Section  4.2  and  Appendix  A.  In  this  section, 
we  investigate  effects  arising  froa  changing  the  ratio  of  inner  spherical 
radius  to  outer  radius. 

Solutions  have  been  obtained  for  an  inner  sphere  of  2 aeters  and  an 
outer  sphere  of  4,  16,  and  22  meters  at.  low  fluences.  Solutions  have  also 
been  obtained  for  an  outer  radius  of  4 and.  10  aeters  at  medium  and  high 
fluences.  The  clear  times  are  marked  on  each  figure,  and  it  is  evident 
when  the  wave  reflected  from  the  outer  sphere  returns  to  the  inner  sphere. 

Figure  4-8  shows  the  surface  current  response  at  90°  on  a 2-meter- 


RT-09939  TIME  (nsec) 

Figure  4-8.  Surface  current  at  90°  cm  a 2-aeter- radius  sphere 
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pulse  width  both  arbitrarily  set  equal  to  20  nsec.  The  snail  outer  sphere, 
with  4-neter  radius,  corresponds  to  a clear  tine  of  only  12  nsec.  The  large 
outer  sphere  corresponds  to  an  84-nsec  clear  time.  Clearly,  the  surface 
current  a^ftlitude  is  enhanced  by  ebout  50%  by  the  interactions  with  the 
snail  outer  sphere.  Note  that  virtually  no  hint  of  oscillatory  behavior 
is  seen. 

Figure  4-9  shows  the  sane  problem  as  that  of  Figure  4-8  except  for  an 
outer  sphere  of  10  neters  and  a higher  photon  fluence,  corresponding  to 


Figure  4-9.  Surface  current  at  90*  on  a 2 -meter- radius  sphere 
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moderate  space-charge  limiting.  (The  peak  emission  current  is  about  4000 
Uftr*}  in  this  enable.)  Dynamic  affects  art  apparent  after  the  occur- 
rence of  each  clear  time.  The  small  outer  ch saber  (R^  • 4 a)  produced 
about  a 50%  amplification  of  the  peak  response,  with  a corresponding  over- 
shoot between  40  and  60  nsec  that  should  significantly  enhance  the  high- 
frequency  spectral  content  of  the  surface  current. 

n 

At  a higher  photon  fluence  level  corresponding  to  4 x 10  tape res, 
the  saae  exaaple  is  repeated  in  Figure  4-10.  In  this  case,  the  charge 
density  just  outside  the  eaitting  surface  is  great  enough  to  produce  plasaa 
oscillations  in  interaction  with  the  rapidly  changing  electric  fields  near 


RT-09941  TIME  (nsec) 

Figure  4-10.  Surface  current  at  90*  on  a 2 -meter- radius  sphere 


the  surface.  !t  appears  that  such  plasma  oscillations  could  btt  producing 
the  shorter  tine  scale  oscillations  in  the  surface  current  curve  in  fig- 
ure 4-10.  The  dynamic  interactions  again  Introduce  about  a $0%  increase 
in  the  peak  surface  current,  and  a substantial  overshoot  between  X$  and 
50  nsec,  in  this  case.  Therefore,  the  medium-  and  hlgh-fluence  dynamic 
interactions  are  qualitatively  the  same  except  for  the  apparent  plasma 
oscillations  in  the  high-fluence  case.  Table  4-1  summarizes  the  dynamic 
effects  shown  in  Figures  4 8,  4-S4 » and  4-10. 


Table  4-1 

SUMMARY  OF  DYNAMIC  INTERACTIONS  BETWEEN  A 2-METER-RAD1US  SPHERE 
AND  A SMALL  Ol/TER  CHAMBER  (R^  « 4 ■) 


Peak 

Surface  Current 

Time  Dependence 

Photon  Fluence 

Amplification 

of  the  Response 

Low 

About  50% 

Monopolar  - no  oscillations 

Medium  (moderate  limiting) 

About  S0» 

Bipolar  - cavity  period 
apparent 

High  (extreme  limiting) 

About  50* 

Bipolar  - cavity  period 
apparent  - plasma  oscilla- 
tions apparently  superim- 
posed on  cavity  ringing 
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S.  COMPARISONS  OF  QUASI-STATIC  AND  FULLY  DYNAMIC  SOLUTIONS  FOR 
ELFCTROMAnNETIC  FIELD  CALCULATIONS  IN  A CYLINDRICAL  CAVITY 

5.1  INTRODUCTION 

The  quasi-static  approximation  has  been  used  extensively  in  the  past 
in  many  calculations  and  computer  programs  for  the  solution  of  electromag- 
netic fields  and  currents  resulting  from  radiation  exposure  of  cavities  of 
different  geometries.  A number  of  qualitative  arguments  have  been  used  to 
justify  the  use  of  the  quasi-static  approximation,  but  no  quantitative 
results  have  been  reported  to  identify  the  range  of  applicability  of  the 
quasi-static  approach  or  the  magnitude  of  the  errors  introduced  by  the 
approximation.  The  fully  dynamic  code  DYNACYL  (Appendix  D)  has  been  applied 
to  a number  of  situations  to  identify  the  ranges  of  validity  for  the  quasi- 
static approximation. 

5.2  DYNAMIC  CODE 

The  dynamic  code  was  developed  tu  treat  a geometry  consisting  of  two 
finite  concentric  cylinders,  each  of  arbitrary  length.  The  fields  are  cal- 
culated for  a cylindrical ly  symmetric  region  enclosed  in  a perfectly  con- 
ducting can.  The  laws  governing  the  generation  of  the  fields  are  Maxwell's 
equations,  and  solutions  are  obtained  by  a finite-differencing  scheme. 

Maxwell's  equations  are  solved  in  cylindrical  coordinates,  assuming 
rotational  symmetry  about  the  cylinder  axis  and  symmetry  as  far  as  direc- 
tion of  rotational  motion  is  concerned.  The  rotational  symmetry  is  similar 
to  assuming  that  all  derivatives  with  respect  to  azimuthal  angle  are  zero. 
However,  this  still  leaves  the  possibility  of  currents  in  the  azimuthal 
direction  (toroidal  currents).  The  assumption  that  there  is  no  preferrejl 
direction  of  rotation  eliminates  this  possibility.  The  fora  of  the  reduced 
equations  is 
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and 
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Yh»*t  equations  expressed  here  in  MKS  units,  ere  solved  for  the  redial 
and  axial  components  of  the  electric  field  (Er,  E^)  and  the  aximuthal  mag- 
netic  field  (M^).  The  medium  in  which  the  fields  are  generated  has  the  same 
permittivity  and  permeability  as  a vacuum.  The  region  of  interest  is  sur- 
rounded by  perfectly  conducting  walls.  It  is  assumed  that  all  fields  are 
initially  sero,  so  the  three  equations  above  imply  the  divergence  equations 

T * Sq  E » p 

and 

V • Uq  H « 0 

if  the  charge  density  P is  determined  by 

* 1 a(rV  Wx 

Jt  ‘ r Jr  “ iz  ’ 


From  the  form  of  the  equations,  it  is  clear  that  they  represent  the 
time  evolution  of  the  fields  driven  by  the  current.  Thus,  it  is  the  cur- 
rent which  must  be  specified  to  determine  the  fields. 


5.3  Q1 JAS 1 -STATIC  APPROXIMATION 

There  are  a number  of  ways  to  view  the  quasi-static  approximation, 
which  ignores  retarded  potentials.  Solutions  are  obtained  by  assuming 
that  the  charge  density  is  quasi-stationary  and  that  the  fields  have 
reached  their  equilibrium  values.  Mathematically,  the  assumption  is 
equivalent  to  the  assumption  that  the  curl  of  the  electric  field  is 
negligib le. 
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Ranges  of  validity  of  tho  quasi- static  approximation  can  bo  estimated 
analytically  by  rearranging  tho  full  tot  of  Maxwell's  aquations  in  terms 
of  tho  oloctric  field  and  source  charts,  nondlmtns ion slicing  tho  resulting 
aquation,  *td  comparing  tons  with  tho  quasi-static  Poisson  aquation.  Two 
important  parsastsrs  smarts  from  such  a coaparison: 

• • v/c 
and 

n • tt/(L/c)  . 

where  v is  the  electron  velocity,  L is  the  cavity  diaension,  c is  the  speed 
of  light,  and  tr  is  the  rise  tine  of  the  e ait ted  pulse  of  electrons.  Thus, 
ft  is  the  ratio  of  electron  speed  to  the  speed  of  li^it,  while  n is  the  ratio 
of  pulse  rise  tint  to  the  tiae  required  for  light  to  traverse  the  cavity. 

fty  requiring  that  the  terms  in  Maxwell's  equations  contributing  to  tha 
fully  dynamic  solution  be  small  relative  to  tha  quasi-static  terms,  inequal- 
ities result: 

n2  » 1 
and 

n >>  6 . 

both  inequalities  must  hold  if  the  problem  is  predominantly  quasi-static, 
because  ft  is  always  less  than  unity,  the  second  inequality  is  automatically 
satisfied  when  the  first  is  satisfied.  The  second  one  is  listed  to  illus- 
trate that  ft  becomes  an  iaportant  parameter  when  the  first  is  not  satisfied. 
Subsequent  results  will  demonstrate  this  behavior. 

These  inequalities  represent  the  conditions  under  which  the  quasi- 
static solution  is  valid;  however,  they  do  not  give  an  indication  of  the 
Magnitude  of  error  involved  when  these  conditions  are  not  met. 

5.4  CODE  RESULTS 

A number  of  fully  dynamic  and  quasi-static  co^iuter  calculations  were 
performed  for  a wide  range  of  parameters  for  which  the  quasi-static  condi- 
tions were  not  always  met.  In  general,  the  dynamic  solutions  vary  from 
the  quasi-static  in  three  important  respects: 
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Peak  amplitude  variations. 

Oscillatory  response  characteristics  of  cavity. 

Delay  time  due  to  finite  velocity  of  light. 

Code  results  were  obtained  for  space-charge- limited  and  non-space-charge- 
limited  cases.  SCL  situations  are  of  great  interest  because  SCL  can  sig- 
nificantly alter  the  time  history  of  the  currents  in  the  cavity  in  a way 
which  cannot  readily  be  ascertained  from  inspection  of  the  emitted  current. 
Thus,  a problem  with  a particular  pulse  time  history  (and  associated  param- 
eter n)  for  which  static  approximation  is  valid  in  the  absence  of  limiting 
may  require  the  fully  dynamic  treatment  when  SCL  occurs.  Calculations  are 
performed  on  a cylindrical  cavity  with  electrons  uniformly  emitted  from  the 
front  face  of  the  cavity,  with  a triangular  time  history  defined  by  the 
full  width  at  half  maximum  (FjVHM)  . Dimensions  and  other  parameters  used  in 
specific  calculations  are  indicated  on  each  of  the  figures.  Mopoenergetio 
electrons  are  used  for  NSCL  calculations , while  Maxwellian  energy  distribu- 
tions are  used  for  SCL  calculations. 

5.4.1  NSCL  Solutions 

5.4. 1.1  Oscillatory  Response.  The  first  set  of  response  calculations  has 
been  selected  to  illustrate  the  tiu>e  history  variations  between  quasi-static 
and  dynamic  calculations  for  several  parameter  values  of  8 and  r,.  Typical 
results  for  the  electric  fields  at  the  front  of  the  cylindrical  cavity  are 
shown  in  Figures  5-1  and  5-2. 

The  conditions  for  the  cylinders  shown  in  Figure  5-1  are  such  that 

2 

8 = 0.27  and  n - 3.  The  inequalities  n >>  1 and  n >>  8 are  well  satisfied, 

and  the  solution  is  essentially  quasi -static. 

The  conditions  for  the  cylinder  shown  in  Figure  5-2  are  such  that 

2 

8 = 0.27  and  n - 0.6.  The  inequelity  n >>  1 is  not  satisfied  and  dynamic 
aspects  of  the  solution  appear.  Note  that  the  dynamic  solution  for  the 
peak  electric  field  at  the  spatial  position  chosen  is  greater  than  the 
quasi-static  solution,  and  that  the  late-time  solution  contains  ringing 
which  is  characteristic  of  the  cavity.  This  trend  is  generally  followed 
by  solutions  for  the  electric  and  magnetic  fields  throughout  the  cavity. 

The  frequency  of  oscillation  is  approximately  equal  to 
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Figure  5-2.  Electric  fields  as  a function  of  time  obtained  from  dynamic  and  quasi-static  codes 

for  parameters  S and  q which  do  not  satisfy  conditions  fot  quasi-static  approximation 
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as  expected. 

A number  of  calculations  were  performed  as  a function  of  the  two  param- 
eters $ and  n to  determine  the  ratio  R of  the  amplitude  of  the  late-time 
oscillatory  response  to  the  peak  amplitude.  This  ratio  is  indicative  of 
the  amount  of  energy  in  the  primary  pulse  which  has  been  coupled  into  the 
cavity.  The  peak  amplitude  of  late-time  oscillatory  response  is  essentially 
constant  in  these  calculations  because  no  damping  has  been  considered.  The 
ratio  R is  shown  in  Figure  5-3  for  a number  of  parameters  0 and  n. 


5.4. 1.2  Response  Delay.  One  shortcoming  of  the  quasi-static  solution  is 
that  retarded  potentials  are  neglected  and,  therefore,  the  response  through" 
out  the  cavity  is  essentially  instantaneous.  This  delay  time  is  not  obvious 
in  the  dynamic  field  solutions  of  Figures  5-1  and  5-2  because  the  fields  are 
calculated  at  the  source  of  the  electron  emission.  The  fields  at  the  back 
of  the  cavity,  however,  are  delayed  by  a tine  approximately  equal  to  L/c, 

as  illustrated  in  Figures  5-4  and  5-5. 

Note  that  when  quasi-static  conditions  prevail,  the  delay  time  is 
insignificant,  as  shown  in  Figure  5-4.  As  the  solutions  become  more 
dynamic  in  nature,  the  delay  becomes  increasingly  more  prominent,  as 
shown  in  Figure  5-5. 

5.4. 1.3  Peak  Amplitude.  In  general,  the  dynamic  solutions  yield  larger 
peak  fields  than  the  quasi-static  solutions.  A parametric  study  of  the 
parameters  0 and  n was  performed  to  determine  deviations  from  the  quasi- 
static solutions.  The  peak  electric  field  at  the  side  of  the  cavity  wall 
was  chosen  as  the  basis  of  comparison.  The  ratio  of  the  peak  electric 
field  obtained  from  the  full  set  of  Maxwell's  equations  is  compared  with 
results  from  the  quasi-static  solution  in  Figure  5-6.  The  trends  are  con- 
sistent with  the  two  inequalities  involving  parameters  0 and  n.  Note  how 
the  results  are  much  more  sensitive  to  0 when  n is  not  much  larger  than 
unity.  Note  also  that  for  small  o,  the  solutions  tend  to  be  quasi-static 
regardless  of  the  electron  velocity.  However,  for  large  n,  the  velocity 
of  the  electron  becomes  important  in  determining  whether  the  solution  is 
quasi-static. 


PEAK  Er  DYNAMIC 
PEAK  £ QUASI -STATIC 


5.4.2  SCL  Solutions 


Dynamic  and  quasi-static  code  calculations  were  also  performed  for 
cases  in  which  the  electron  currents  were  space-charge  limited.  The  fol- 
lowing calculations  are  intended  to  demonstrate  how  a situation  with  a 
particular  pulse  shape  can  be  essentially  quasi-static  when  there  is  no 
SCL  but  may  require  a full  dynamic  treatment  when  SCL  occurs.  This 
situation  arises  because  the  fields  alter  the  electron  trajectories  in 
such  a way  that  the  effective  current  pulse  shape  at  locations  within 
the  cavity  is  significantly  different  from  the  emission  current  pulse 
shape.  The  frequency  content  associated  with  the  emission  current  pulse 
shape  might  well  satisfy  the  quasi-static  requirements  for  8 and  n,  but 
:he  effective  current  pulse  shape  which  determines  cavity  fields  may  not. 
Such  a case  is  treated  in  this  section. 

Electrons  are  emitted  uniformly  from  one  face  of  a pillbox  cavity  with 

a length-to-diameter  ratio  of  unity.  The  emitted  electron  time  history  is 

triangular,  with  rise  and  fall  times  equal  to  S nsec.  The  emitted  electron 

energy  distribution  is  proportional  to  exp{- | E-Eq | /EQ)  where  E^  is  selected 

such  that  the  average  electron  energy  is  0. 1 MeV.  Emission  current  levels 

-» 

ot  0.1,  1.0,  and  10  amp/cm"  were  chosen,  roughly  corresoonding  to  small, 
medium,  and  large  a.-tounts  of  SCL. 

The  amount  of  SCL  can  be  quantitatively  gauged  by  the  fractional  amount 
of  emitted  current  which  reaches  the  opposite  side  of  the  cavity,  as  shown 
in  Figure  5-7  for  each  of  the  three  cases. 

Based  on  the  emitted  current,  the  parameters  n2  and  8 have  values  of 
9 and  0.55,  and  the  solutions  might,  therefore,  appear  to  be  essentially 
quasi-static.  The  electric  fields  at  the  back  of  the  cavity,  as  shown  in 
Figures  5-8,  5-9,  and  5-10  for  the  three  cases,  indicate  that  the  solution 
is  quasi-static  only  for  the  NSCL  solution,  and  that  the  dynamic  aspects 
of  the  solution  become  more  apparent  as  the  amount  of  SCL  increases.  This 
result  is  not  surprising  since  the  fields  throughout  the  majority  of  the 
cavity  are  more  dependent  on  the  transmitted  current  than  on  the  emitted 
current.  The  effect  of  SCL  is  to  decrease  the  effective  pulse  rise  time. 

It  would,  therefore,  appear  that  an  effective  value  of  n should  be  defined 
based  on  the  rise  time  of  the  transmitted  current  rather  than  on  the  emitted 
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Figure  5-9.  Electric  field  at.  back  versus  lim  for  Medium  SCL  case 


current.  Effective  values  of  8 and  n based  on  the  transacted  current  are 
shown  in  Table  S* 1 . It  is  evident  that  the  solutions  should  be  essentially 
quasi ^static  for  the  first  two  cases  but  not  for  the  third. 


Table  5-1 

PARAMETERS  8 AND  n FROM  RESULTS  OF  THREE  DIFFERENT 
EMISSION  CURRENT  RUNS  WITH  SPACE -QIARGE- LIMITING 


Emission 
Current 
Level 
(amp/cnr 1 

8 

R<  se 
Time 

(nsec) 

nEFF 

n*EFF 

nBFF/a 

0.  1 

0.55 

5.0 

3.0 

9.0 

16 

1.0 

0.55 

2.8 

1.6 

2.6 

4.7 

10.0 

0.55 

1.4 

0.82 

0.66 

1.2 

8 is  the  average  value  of  electron  velocity  over  the  speed  of  light. 
n£FF  *s  r*se  time  of  transmitted  current  pulses  divided  by  L/c  time  cf  the 

cavity. 

Quasi-static  solutions  require  o2  '*>  1;  n/8  >>  1. 

5.5  SUMMARY  AND  CONCLUSIONS 

A fully  dynamic,  self-consistent  IEM*  code,  DYNACYL,  has  been  developed 
and  exercised  to  determine  the  range  of  validity  of  previously  developed 
quasi-static  calculations.  Solutions  for  the  dynamic  treatment  vary  from 
the  quasi-static  approach  in  three  important  respects.  The  dynamic  solu- 
tions generally  display  higher  peak  field  values  than  quasi-static  solutions, 
are  delayed  in  time  due  to  the  finite  speed  of  light,  and  contain  oscilla- 
tory late-time  solutions  characteristic  of  resonances  associated  with  the 
particular  geometry  of  interest. 

Two  parameters  have  been  defined  which  determine  whether  the  solution 
to  a particular  problem  will  display  dynamic  effects  or  will  be  essentially 
quasi-static.  It  has  been  found  that  problems  which  are  essentially  quasi- 
static without  SCL  may  display  highly  dynamic  behavior  under  SCL  conditions. 
The  SCL  causes  fields  and  currents  to  reach  maximum  values  much  sooner  than 
the  emitted  pulse  and  thereby  increases  the  frequency  content  of  the  effec- 
tive driving  function. 
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Finally,  it  has  been  determined  that  the  quasi-static  solutions  yield 
essentially  the  same  results  (to  within  a factor  of  two)  as  dynamic  solutions 
over  virtually  the  entire  range  of  interest  of  IEMP  problems,  with  the  not- 
able exception  of  the  late-time  response  after  the  end  of  the  emitted  pulse. 
This  portion  of  the  response  is  electromagnetic  in  nature  and  cannot  be 
determined  from  the  quasi-static  approach. 
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APPENDIX  A 

SCALING  SGEMP  EXPEKlfCNTS* 
1.  INTRODUCTION 


! 


In  many  system-related  experiments  dealing  with  SGEMP  effects,  it 
may  be  necessary  or  desirable  to  scale  one  or  more  of  the  many  variables 
available  to  the  experimenter.  Insight  into  scaling  laws  is  necessary  to 
design  meaningful  experiments,  choose  proper  excitation  sources,  and  relate 
results  to  systems  of  interest.  The  effort  described  in  this  n-«te  is  a 
first  attempt  at  putting  down  a number  of  obvious  relationships  between 
scaling  parameters. 

We  begin  first  by  listing  the  important  quantities  which  can  most 
readily  be  varied  by  the  experimenter: 

t = time  (varied  by  changing  the  time  history  of  the  photon  pulse) 
r = dimension  (characteristic  dimension  of  test  object) 

Y = incident  photon  fluence 
v = electron  velocity 

There  are  also  a number  of  ancillary  quantities  which  cannot  read- 
ily be  varied: 

q = electron  charge 
m = electron  mass 
e = permittivity  of  free  space 
p = permeability  of  free  space 
a - conductivity 

We  seek  information  on  how  the  important  variables  of  time,  dimen~ 
sion,  fluence,  and  electron  velocity  will  affect  the  electromagnetic  quanti- 
ties of  interest,  which  include: 

E = electric  field 
H = magnetic  field  strength 
1 = current 
p = charge  density 
J = current  density 
AV  = potential  difference 
*Work  performed  under  I RT- sponsored  IR&D  program. 


Pmife!  pip  kM 


i'm 
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It  is  convenient  to  relate  the  scaled  quantities  denoted  by  a (') 
from  the  unsealed  quantities  with  the  following  relationships: 
t'  « t/T 
r'  « r/R 
y#  » rx 

v'  = 8v 

where  the  factors  r,  R,  r,  and  8 are  the  dimensionless  scaling  factors.  We 
assume  that  q,  m,  u,  anJ  c do  not  scale  and  that  a is  infinite. 

In  this  note,  we  determine  the  scaling  laws  for  the  linear  prob- 
lem (Section  2)  and  the  nonlinear  problem  (Section  3).  Our  solutions, 
strictly  speaking,  apply  to  the  space  surrounding  perfect  conductors,  as 
well  as  to  currents  flowing  on  the  conductors.  Coupling  to  cables  and 
objects  with  dielectrics  is  briefly  discussed  in  Section  4.  Solutions  for 
several  example  problems  using  the  nonlinear  DYNASPHERE  SGEMP  code  are  pre- 
sented in  Section  5 for  illustrative  purposes.  Finally,  system  implica- 
tions are  discussed  in  Section  6,  the  summary. 


2.  LINEAR  SCALING 


In  this  section  we  apply  the  scaling  laws  to  Maxwell's  equations 
to  determine  how  the  electromagnetic  quantities  of  interest  scale  for  linear 
problems  in  which  the  electron  motion  is  not  perturbed  by  the  fields.  The 
emitted  current  becomes 


j'  _ r dy  p dy 

Je  - f ~3C  ~ f Ft  dt  ’ 


(1) 


The  charge  density  at  the  emitting  surfaces  located  at  point  Rq  for  each 
velocity  increment  becomes 


J' 

e 


p'(R<J)  v'  , 


p;cr<P 


J - 
e 

V' 
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The  scaled  equation  of  notion  for  the  electron  in  the  linear 
region,  where  the  electron  trajectory  is  unaffected  by  the  electron  and 
magnetic  fields,  is 

r'(t ')  * v't'  ♦ Rq'  , (2A) 

where  R'  is  the  position  of  the  electron  on  the  emitting  surface.  Substitut- 
ing for  v',  t',  and  Rq, 

r-(f)  - . 

The  solution  scales  if  the  following  relationships  hold. 

j - R (3) 


re  ■ S tyt  * V * TT  • (4) 

The  position  of  the  electron  scales;  therefore,  the  ratio  of 
the  charge  density  at  a point  in  space  to  the  charged  density  at  the  emit- 
ting surface  in  the  scaled  geometry  equals  the  same  ratio  in  the  unsealed 
geometry  for  electrons  of  velocity  v. 

W Pv(R0} 

Pv(r')  * py(r) 

From  Eq.  2, 

Pyffyp  rT 

■ »v(r)  TTrJ-  * T ■ (5) 

This  result  holds  for  charge  as  a result  of  any  emitted  velocity 
increment  and,  therefore,  holds  for  all  velocity  increments.  The  total 
charge  density  due  to  all  velocity  increments  scales  properly  provided  all 
velocity  increments  scale  as  g (the  velocity  or  energy  distribution  remains 
unchanged).  ~hus, 

p'(r-)  = p(r)  y-  . (6) 
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IT  -1  dll 
n*^‘T3t 


Therefore, 


R2E* 

Pt2 


E , 


or 


E'.^G  . 


Using  condition  3,  3 = t/R, 


E " = 


TtSE 

R 


(12) 


The  expression  for  the  scaled  electron  field  given  by  Eq.  12  is 
inconsistent  with  Eq.  8 unless 


8 = 1 . 


(13) 


Thus,  we  see  that  proper  scaling  can  be  achieved  only  if  the  electron  velo- 
city is  not  scaled. 

The  potential  difference  between  any  two  points  is  given  by 


fy-r' 
'Cr')  = - / E- 

y - *n 


AV 


(y)  Jy 


(14) 


J 


r y 3 r/R 


E '(y)  dy 


y » Ro/R 


ry  ° R 

> = R, 


E'(y)  ^ 


o 

y = R 


t/ 

' - Ro 


¥ 7 


AV'(r')  = ~ AV(r) 
R 


(15) 


'V 


V 
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The  total  currents  scale  as  follows. 


Condition  13,  along  with  condition  3,  requires 


(16) 

(17) 


T * R £ a , (18) 

where  a is  defined  as  the  primary  scaling  factor  for  both  time  and  space. 
Scaling  results  are  summarized  in  Table  B-l . 


Scaling  Factors 

t*  * t/t 
r‘  * r/R 
v'  = 8v 
y'  * rY 


Table  A-l 

Scaling  Requirements 

#4  R * t 
8 = 1 


Scaling  Results 

j'  = or  j 

e 

p'  = arp 
E ' * TE 

i r = m 

AV'  = - AV 
a 


In  summary,  scaling  is  possible  for  a linear  matrix  of  electrons 
in  free  space  in  the  presence  of  perfect  conditions  if  the  relationships  in 
Table  A-l  are  satisfied.  The  varioi s quantities  will  scale  as  shown.  The 
development  here  does  not  treat  dielectrics  or  imperfect  conductors. 

If  we  require  that  the  potential  at  the  scaled  point  r'  be  the 
same  as  that  at  the  unsealed  point  r,  then  we  require  that 


r * a , 


(19) 


a condition  which  vc  will  show  in  Section  3 is  necessary  for  proper  seeling 
of  the  nonlinear  solution. 


3.  NONLINEAR  SCALING 

Scaling  for  cases  in  which  the  electron  trajectories  are  perturbed 
by  the  fields  can  be  treated  by  replacing  Eq.  2A  by 

rC»t' 

I F[r'(0,  C]  dCdn  v'(r',«  dt  ♦ R^  (20) 

n*0  ^C«0  *'c*0 


where 


F'  « q(E  ' + v'  B')  . 

We  first  replace  t'  by  t/x  on  the  right  side: 

J.n-t/x  Can  rt^t/T 

F[r'(0»€]  dEdn  ♦ I v'(r'.C)  dK  + Rq 
n=0  ^*0  *'C«0 

Now  replace  the  arbitrary  integration  variable  n by  n/x: 


Jr* n=t  pK-T\/T  /*t= t/T 

I F'[r'(0,ej  dC^Uj  v-(r-.C)  dc  ♦ . 

n=0  *'C»0  ,'C*0 

Next  replace  the  arbitrary  integration  variable  £ by  i/x  and  the  coordinate 
R0  b*  R0/R: 


C=t/r 


OS-n 

>o 


F'(r'(C/n).e/n] 


Using  the  expression  for  E'  and  H'  in  Table  A- 1,  Eq.  21  becomes 
F*  * q(E'  ♦ v'  X 8')  = q(TE  + Tv  x B)  » TF  , 
and  the  equation  of  motion  becomes 


r'(t)  » 4 f f F d*  dn  * 

T *0  •'o 

The  position  r'(t)  Will  scale  properly  if 

X 1 . I 

2 * t R ' 

T 

But  t ■ R « o,  so  that 

r „ 1 . ? 

2 a a 

01 


Therefore,  for  the  position  of  the  electron  to  scale  properly,  the  follow- 
ing relation  must  hold. 


1 F • 

This  condition  is  equivalent  to  condition  19,  requiring  that  the 
potential  at  the  scaled  point  r'  equal  that  at  the  unsealed  point  r. 

We  summarize  the  scaling  requirements  for  the  linear  and  nonlinear 
cases  in  Table  A-2. 


Table  A-2 


Scaling  Parameters 

Scaling  Factors 

t ■*  » t/r 

x * a 

r " = r/R 

R = a 

v'  = fiv 

B * 1 

Y'  = Ty 

r = a 

Scaled  Field  Quantities 

.1'  = a2  J 
2 

P * a p 
E'  = aE 
H'  = aH 
Av'  = Av 
I " = I 
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4>  SAMPLE  PROBLEMS 


To  verify  the  scaling  laws  discussed  in  the  previous  several 
chapters,  we  have  exercised  the  DYNASPHERE  SGEMP  code  for  several  problem 
of  interest.  The  conditions  for  the  two  sets  of  problems  shewn  in  Table 
A>3  were  chosen  so  that  one  set  would  be  linear  and  one  set  would  be  spnee- 
chargc- limited  and.  therefore,  nonlinear. 

Magnetic  fields,  electric  fields,  and  surface  currents  on  a spher 
ical  surface  were  calculated  at  the  following  points. 

E St  6 » O' 
r 

B at  6 * 90* 

♦ 

I through  a plane  at  Q * 90® 

Results  are  shown  in  Table  A-4  and  graphically  illustrated  in  Figures  A-l 
through  A- 6. 

Note  that  the  results  are  plotted  as  a function  of  E x R,  H x R, 
and  I versus  t/(2*R/c)  for  a triangular  pulse  of  rise  time  and  full 
width  at  half  maximum  of  irR/c.  While  the  scales  at  first  seem  abstract, 
the  results  are  quite  general  in  that  they  apply  to  any  sphere  with  radius 
R excited  by  a pulse  with  rise  time  xR/c. 


Table  A-4 


TABULATED  RESULTS  FOR  TWO  PROBLEM  SETS 


Time 

(nsec) 

J 

e 

(a«p/m2) 

E (0*) 
(volt/m) 

H (90*) 
(amp/m) 

Surface 
Current  (90*) 

(a*p) 

Inner  radius  0.2  a; 

peak  emission  current 

• 10  aap/a2;  t 

• 2 nsec 

• 2 r 

• 2 

1 

4.7S 

2.736  x 10* 

-6.006  x 10 

-7.454  x 10 

2 

9.7S 

1.063  x 103 

-2.668  x 10"1 

-3.3112  x 10*1 

3 

S.2S 

1.714  x 103 

-4.S6S  x 10" 1 

-5.6660  x 10* 1 

4 

2.S  x 10*1 

1.646  x 103 

-4.643  x 10'1 

-5.7627  x 10“l 

S 

0 

1.271  x 103 

-3.349  x 10'1 

-4.1566  x JO'1 

Inner 

radius  ■ 2 «; 

peak  emission  current 

0.1  amp/m2;  tr  • 

20  nsec 

10 

4.7S  x 10‘2 

2.736  x 101 

-6.004  x 10"3 

-7.4522  x 10'2 

20 

1 x 10‘1 

1.066  x 102 

-2.670  x 10'2 

-3.3140  x 10'1 

JO 

S.25U0  x 10' 

-2 

1.716  x 102 

-4.573  x 10‘2 

-5.6753  x 10’1 

40 

2.S  x io~3 

1.648  x 102 

-4.648  x IO*2 

-5.7689  x 10'1 

50 

0 

1,273  x 102 

-3.357  x 10‘2 

-4.1666  x IO'1 

4 2 

Inner  radius  « 0.2  m;  peak  emission  current  • 5 x 10  amp/»  ; tr  ■ 2 nsec 


4 

5 

1 

1 

0.5 

1.1875 

X 

10 

3. 257 

X 

10 

-4.962  x 

io1 

>6.584  x 10 

1 

2.4375 

X 

io4 

1.291 

X 

io6 

-2,335  x 

io2 

-2.8982  x 10 

1.5 

3.6875 

X 

io4 

2.019 

X 

io6 

-2.672  x 

102 

-3.3169  x 10' 

2 

4.9375 

X 

io4 

1.961 

X 

io6 

-1.502  x 

io2 

-1.8642  x 10 

2.S 

3.8125 

X 

io4 

1.838 

X 

io6 

-1.114  x 

io2 

-1.3831  x 10 

3 

2. 5625 

X 

io4 

1.625 

X 

io6 

-8.582  x 

IO1 

-1.0652  X 10 

3.5 

1.3125 

X 

io4 

1.306 

X 

io6 

-1.394  x 

io2 

-1.7301  x 10 

4 

6.25  x 

102 

7.963 

X 

10s 

-7.442  x 

io1 

-9.2373  x 10 

4.5 

0 

4.54  x 1 

LO5 

1.251  x 

io2 

1.5523  x 102 

5.0 

0 

3.937 

X 

10S 

1.950  x 

io2 

2.4201  x IO1 
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Table  A-4  (continued) 


* 


t- 


TiM 

(nsec) 

J 

e 

(amp/m2) 

E (0*) 
(volt/m) 

M (90*) 
(amp/m) 

Surface 
Current  (90*) 

(amp) 

t Inner 

radius  * 2 i; 

peak  emission  current 

2 2 
■ 5 x 10  amp/m 

; tr  * 20  nsec 

S 

1.87S  * 102 

3.25?  x 104 

-4.9561 

-6.5464  x 10*1 

10 

2.437S  * 102 

1.291  x 105 

-2.334  x 10l 

-2.8972  x 10 

IS 

3.6875  * 102 

2.019  x 105 

-2.672  x 1C1 

-3  3168  x 102 

20 

S x 102 

1.963  x 10S 

-1.502  x 101 

-1,8645  x 102 

2S 

3. 8125  x 102 

1,838  x 105 

-1.118  x 101 

-1.3882  x 102 

50 

2. So2S  x 102 

1.626  x 10S 

-8.492 

-1.0540  x 102 

3S 

1.3125  x 102 

1.306  x 105 

1.401  x 101 

-1.7393  x 102 

40 

6 . 25 

7.962  x 104 

-7.441 

-9.235  x 10'1 

4S 

0 

4. 54  >;  104 

1.261  x 101 

1.5655  x 102 

SO 

0 

3.93S  x 104 

1.933  x 101 

2.3993  x 102 

I 

! 

i 

t 

| 

i 


Figure  A-2  Radial  electric  field  at  the  spherical  surface  (0  = 90°)  for  the 
two  low-emission  cases.  Electron  emission  is  sufficiently  low 
that  the  problem  is  linear.  Excitation  pulse  rise  time  tr 
is  irR/c. 
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Figure  A-5  Radial  electric  field  at  the  spherical  surface  (0  * 90°)  for  the 
two  high-emission  cases.  Electron  emission  is  sufficiently  high 
that  the  problem  is  nonlinear.  Excitation  pulse  rise  time  tj  is 


2 meters 


Figure  A-6  Surface  emission  passing  through  a plane  located  at  0 * 90°  for  the 

two  high-emission  cases.  Electron  amission  is  sufficiently  high  that 
the  problem  is  nonlinear.  Excitation  pulse  rise  time  tr  is  wR/c. 


S.  SUMMARY 


Of  the  many  variables  available  to  the  SGENP  experimenter,  the 
ones  of  most  interest  arc  the  object  dimension,  the  pulse  time  history,  the 
incident  fluence.  and  the  incident  photon  energy  spectrum  (or  t>ie  emitted 
electron  velocity  distribution).  In  this  note  we  have  explored  the  possi- 
bility of  changing  one  or  more  of  these  parameters  to  values  more  conven- 
ient to  the  experimenter,  while  retaining  the  essential  features  and  results 
corresponding  to  the  original  parameter  values. 

First,  we  find  that  one  cannot  change  the  electron  velocity  and 
hope  to  scale  the  solutions.  A simple  "Gedanken"  experiment  illustrates 
the  problem.  Consider  a pillbox  with  steady  current  flowing  from  one  face 
to  the  other.  The  magnetic  field  within  the  cavity  is  proportional  to  the 
emitted  current  density  and  is  independent  of  the  velocity.  On  the  other 
hand,  the  electric  field  is  proportional  to  the  charge  density  which,  in 
turn,  is  inversely  proportional  to  the  electron  velocity.  Thus,  the  mag- 
netic and  electric  field  ratios  change  as  a function  of  electron  velocity, 
and  the  solutions  cannot  scale. 

Next,  we  find  that  the  dimensions  and  pulse  time  must  scale  by 
the  same  factor.  This  scaling  is  required  to  keep  the  ratio  of  excitation 
wavelength  to  object  wavelength  constant.  It  is  also  required  to  ensure 
that  the  electron  reaches  the  scaled  distance  from  the  object  at  the  cor- 
rect scaled  time. 

Finally,  we  find  the  total  charge  emitted  must  be  changed  by  the 
inverse  factor  used  for  pulse  and  object  dimension  if  we  are  to  retain  the 
same  potential  at  the  scaled  coordinate.  For  nonlinear  problems  in  which 
space-charge  limiting  is  important,  it  is  essential  to  have  the  potential 
difference  between  two  scaled  points  equal  to  that  between  two  unsealed 
points  (along  the  same  path). 

We  therefore  conclude  that  the  essential  features  of  a problem 
are  unchanged  if 

1.  the  electron  velocity  is  constant, 

2.  the  pulse  time  history  is  scaled  by  a. 
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3.  the  object  dimension  is  scaled  by  a, 

4.  the  incident  fluence  is  scaled  by  1/a. 

Such  a scaling  will  result  in  the  following  changes  in  field  values,  assum- 
ing electron  motion  in  free  space  in  the  presence  of  perfectly  conducting 
boundaries. 


t*  * t/a 
r'  * r/a 
r'  * ar 
v'  * v 

t*  2 . 

J * a J 


, 2 
P 8 a p 

E ' * aE 

H'  = aH 

AV'  » AV 

I'  = I 


By  scaling  time,  we  mean  that  ail  features  of  the  incident  photon 
pulse/emitted  electrons  change  by  the  factor  a:  the  pulse  is  essentially 

expanded  or  compressed  by  the  factor  a. 

Scaling  problems  associated  with  the  electron  velocity  are  of  par- 
ticular concern  because  of  the  wide  variety  of  photon  sources  currently 
available  or  under  consideration.  The  results  of  this  note  indicate  that, 
strictly  speaking,  one  cannot  modify  other  variables  such  as  object  dimen- 
sions, machine  output,  time  history,  etc.,  to  offset  changes  in  electron 
velocity.  Therefore,  one  cannot  directly  extrapolate  results  from  one 
photon  energy  range  to  another.  However,  we  note  that  the  electron  velo- 
city has  a square-root  dependence  on  energy,  and  therefore,  the  results 
are  not  extremely  sensitive  to  changes  in  the  energy  distribution. 

It  should  be  noted  that  decreasing  an  object  size  by  the  factor 
a requires  a corresponding  decrease  in  the  pulse  time  history  by  the  factor 
a;  however,  the  photon  fluence  need  not  be  increased  by  a corresponding 
factor  a if  the  problem  is  in  the  linear  regime  where  space-charge  limiting 
is  unimportant. 
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APPENDIX  B 
DYNASPHERE 


This  appendix  contains  a mathematical  description  of  the  code  .)YNA- 
SPHERE,  which  has  been  developed  to  numerically  evaluate  electron  motion 
and  electromagnetic  field  generation  in  a region  bounded  by  two  concentric 
spheres.  The  spheres  are  taken  to  be  perfectly  conducting,  and  the  prob- 
lem characteristics  are  taken  to  be  rotationally  symmetric,  reducing  the 
problem  to  two  dimensions.  The  electromagnetic  fields  are  calculated  from 
the  full  set  of  Maxwell's  equations  in  this  simplified  geometry,  and  the 
electric  fields  are  used  to  influence  the  electron  motion. 

Electron  emission  from  the  inner  sphere  is  the  source  term  for  driv- 
ing the  problem.  Currents  and  fields  between  the  two  spheres  are  the  quan 
tities  which  result  from  the  calculation.  The  electron  emission  must  be 
specified  in  space  and  time.  For  example,  in  the  case  where  the  emission 
is  due  to  photon  interaction  with  materials,  the  photon  energy  and  time 
spectrum  determines  the  emission  characteristics  of  the  electrons.  (Par- 
ticles are  used  to  represent  large  numbers  of  electrons.)  The  quantities 
calculated  directly  include  the  currents  in  the  region  between  the  boun- 
daries as  well  as  the  electric  and  magnetic  fields  in  that  region.  Apply- 
ing Maxwell's  equations  properly  at  the  boundaries  gives  surface  currents 
and  charge  densities. 

The  remainder  of  this  appendix  outlines  the  method  of  solving  the 
field  equations,  the  conversion  of  particle  motion  int^  currents  useful 
for  driving  these  equations,  and  the  method  of  calculating  the  particle 
motion.  Since  the  first  two  of  these  procedures  is  given  for  a general- 
ized coordinate  system,  their  description  is  given  initially.  The  spe- 
cialization to  the  particular  coordinate  scheme  used  here  is  then  dis- 
cussed, and  the  particle  motion  is  treated  in  these  coordinates. 

The  numerical  solution  of  Maxwell's  equations  is  straightforward  in 
a region  enclosed  by  perfectly  conducting  surfaces,  such  as  the  region 
considered  in  DYNASPHERE.  The  initial  conditions  are  that  all  fields 
and  charge  densities  are  zero.  The  currents,  as  calculated  from  the 
motion  of  the  charged  particles  injected  into  the  region  (see  below). 


»re  the  quantities  which  drive  the  time  evolution  of  the  fields.  Thus, 
for  a medium  with  the  permittivity  (cQ)  and  permeability  (tig)  of  free 
space,  Maxwell's  equations  reduce  to 
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with  the  initial  conditions 

t « 0 at  t » 0 , 

H » 0 at  t ■ 0 , 

E x n ■ 0 where  n is  the  normal  to  the  bounding  surface  . 

Putting  these  equations  into  numerical  form  for  solution  is  also 
straightforward,  and  may  be  carried  out  directly  in  two  dimensions.  How- 
ever, by  performing  the  task  in  three  dimensions,  one  is  forced  naturally 
into  a symmetry  which  is  very  convenient  and  not  obvious  in  the  two- 
dimensional  case. 

To  complete  the  geometrical  generality  (and  permit  a si~ple  method 
of  varying  zone  spacing! , the  space  under  consideration  is  taken  to  be 
metrized  by  the  generalized  orthogonal  coordinates  (qj,  q?,  q3)  where  the 
order  is  such  that  the  coordinate  system  is  right-handed.  Using  the  nota 
tion  of  Margenau  and  Murphy,  displacements  in  real  space,  -s,  may  be 
related  to  displacements  in  q space  by  the  functions  Q,  where 

ds  = dq.  . 


Using  this  definition  and  the  definition  of  the  curl,  Vx,  as  the  path 
integral  of  a quantity  in  the  left-hand  direction  (Cauchy  rule)  around 
a closed  path  divided  by  the  area  of  the  enclosed  surface.  Maxwell's 
equations  in  finite-difference  form  can  be  represented  by 


dr 

H.  Margenau  and  G.  M.  Murphy,  "The  Mathematics  of  Physics  and  Chem- 
istry," Princeton,  D.  van  Nostrand  Company,  Inc.  (1956). 
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R,  (i+'i.j.k)  * H.  (i*‘a,j,k)  - (i+’a.j.k) 

* i 0 

J . At^J I 

+ 1C0  Q2  A^2  Q3  j >k)  Aq^[ 


x J^‘5  (i*S.M.k)  Q3  (i+,i» j+,i,k)  Aq5  - (i^aj-'a.k)  Qs  ( Aq 

- [*»'>  (i+'j.j  ,k+'i)  Q,  (i+’j.j  .k-*^)  Aq,  - H2+i  (i*'a,  j ,k-5i)  Q2  (i+’a,  j ,k«*j)  Aq 


VJ+V2  (i.j+'a.k^a)  = h"*'j  (ij  + 'a.k*',) 


VQ  Q-,  (i.j**a.k+‘i)  Aq,  Q3  ( i J + 'a.k+'j)  Aq3 


X|[E5+1  ti.j  + ltk*‘a)  Q3  (i  ,j  + l ,k  + 'j)  Aq3  - 1 (i.j.k+’a)  Q3  (L.j.k  + 'j)  Aq.J 

- je2+1  (i.j*‘a,k+l)  Q2  (ij+'a.k+l)  Aq,  - r"+ 1 (i.j+'j.k)  Q,  (i,j  + l?,k)  Aq,JJ 

with  the  continuity  equation  (not  essential  to  determining  the  fields,  hut 
of  interest  for  itself  and  for  converting  particle  motion  into  equivalent 
currents)  represented  by 

P (i,J ,k)  = p (i ,j ,k) 


Qj  ( i , j ,k)  Aqj  Q2  (i,.i,k)  Aq2  Q3  (ij.k)  Aq. 


x (i+Js*j»k)  Q,  (i+xi, j ,k)  Aq2  Q3  (i+'i.j.k)  Aq3 

- (i-’a.j.k)  Q2  ,k)  Aq2  Q-  (i-’ij.k)  AqJ 

♦ ^2+  J (i»j+12»k)  Qj  (i,.i+li,k)  Aq1  Q3  (i  i , j+li,k)  Aq. 


- (i  »i“-3*k)  0,  (ij-'a.k)  Aq 


, Q,  (i»j-'i.k)  Aq-l|  . . . 


♦ 


Jj4*  (i  * j Qj  Aq1  Q2  Aq2 

» 

“ Jj4j  (i,j,k-V)  Qt  (i.j,k-*i)  &ql  Q2  (i,j,k-*s)  Aq, 

The  notation  used  in  these  equations  requires  some  comment.  The  super- 
script refers  to  the  tine  step.  Thus,  some  quantities  are  centered  in  time 
and  some  are  at  boundaries  in  time.  Corresponding  to  this,  there  are  two 
time  steps:  the  time  step  connecting  quantities  centered  in  time,  Atn , 

and  the  time  step  connecting  quantities  at  boundaries  in  time,  At  \ 

Subscripts  refer  to  directions  in  the  generalized  coordinate  space. 
Thus,  is  the  component  of  the  electric  field  along  the  direction  of  a 
displacenwnt  in  space  given  by  a displacement  in  qj  at  the  spatial  point 
in  question.  The  quantities  in  parentheses  refer  to  the  position  in  space. 
Thus,  Ej(a,0,Y)  is  evaluated  at  the  point  in  space  determined  by  the 
coordinates 

ql  “ a Aql  * q 1 (MIN)  * 
q2  = 3 Aq2  ♦ <I:(MIN) 

q3  = Y Aq3  * q3(MIN  * 

where  the  minimum  value  of  the  coordinate  is  specified  for  convenience, 
allowing  the  spatial  boundaries  to  be  other  than  zero  in  the  q space. 

It  will  he  noted  that  the  grid  spacing  in  q space  is  uniform.  Fur- 
ther, if  one  of  the  Q's  is  zero  at  a point  of  interest,  the  procedure 
fails.  In  fact,  at  such  points  the  coordinate  system  does  not  met  rise 
real  space.  The  failure  is  that  many  points  in  q space  correspond  to 
one  point  in  real  spate.  Such  cases  must  be  treated  specially. 

The  equations  for  ihc  other  components  of  the  electric  and  magnetic 
fields  are  obtained  by  cyclically  permuting  the  integer  subscripts  and 
the  corresponding  coordinates  a,  is,  and  y. 

Description  of  Charged-Particle  Motion  by  Currents 

In  the  numerical  solution  of  Maxwell's  equations,  currents  are  eval- 
uated at  discrete  points  in  space  and  time.  Particles,  representing 
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electrons,  move  in  a generalized  coordinate  system  in  time.  An  interpola- 
tion scheme  must  be  constructed  to  translate  the  continuous  motion  of  the 
particles  into  a set  of  discrete  currents.  The  scheme  chosen  is  such  that 
the  residual  charge,  as  calculated  from  the  time  integral  of  the  divergence 
of  the  current,  is  2ero  after  a particle  has  passed  into  and  then  out  of 
a region.  (This  does  not  hold  for  boundary  zones  since  currents  outside 
the  region  of  interest  are  not  considered.) 

In  the  present  zoning,  currents  are  evaluated  at  zone  centers  in  their 
own  direction  and  zone  boundaries  in  perpendicular  directions.  The  current 
is  time-centered.  Accordingly,  as  forced  by  the  continuity  equation,  the 
charge  densities  are  at  zone  boundaries  in  space  and  time. 

The  present  objective  is,  then,  to  take  a particle  from  point 


a a a .a. 


to  the  point 

, b b b h 
(qrq2,q3,t  j 

where  t = t + At . That  is,  the  currents  representing  such  h translation 
must  be  generated. 

Since  the  grid  spacing  is  uniform  in  the  q's,  linear  interpolation 
can  be  performed  in  the  transverse  coordinates.  This  can  be  seen  from 
the  formulation 

a At  K-q  I 1 

dj  (i  + b,  j ,k  ,t J ► ^-)  = Qp  I Q i (i  + 'i.j.k) 

[Qj  (i*\s,j,k)  Aqj  Q2  (i  + 'j.j.k)  Aq-,  Q5  (i+'j.j.k)  AqJ 
q2  " q2(-j)  + (q2  " q2^  ^ ' t /At  |”| 

q3  ’ q3^  + ^q3  ” - t J)  /At 

Aq3 

for  point  (q^q^.qj)  = [(i+’.jAqj , jAq^;  .kAq^]  , particle  of  charge  Qp.  Mere, 
is  the  current  in  the  q^  direction. 


S3 


* > 


The  first  bracket  contains  the  real-space  particle  velocity  in  the  qj 
direction.  The  second  bracket  contains  the  real-space  volume  element . The 
integral  represents  the  time- average  fractional  distance  from  the  point  of 
interest  in  the  transverse  directions.  It  can  be  seen  that  the  particle  is 
treated  as  a volume  element  in  q space  of  dimensions  Aq } , fiq, , Aq^,  and  the 
integral  represents  the  transverse  area  overlapping  the  zone  of  interest. 

The  time  limits  are  such  that  they  are  within  the  range  of  interest 

k »a 

t , t * and  the  particle  is  within  one  zone  of  the  point  of  interest  in 
its  linear  traversal: 

f a . b a.  (t  - ta)  1 
Ml  » [q4  * 

q,  * , . 

;u  * ... 

a 

from  point  a to  point  b. 

In  practice,  a particle  is  tracked  from  point  a to  point  b in  inter- 
vals of  zone  crossings  so  that  the  appropriate  values  of  i,  i,  and  k are 

easily  established.  The  treatment  for  .1,  and  .J,  is  related  to  that  of  J, 

2 a 1 

by  the  same  symmetry  as  mentioned  in  the  discussion  of  the  field  calculation. 

Restriction  to  Two-Dimensional  Spherical  Coordinates 

Here  the  application  to  the  special  case,  where  the  coordinates  are 
radius  r and  polar  angle  0,  is  considered.  Azimuthal  velocities  of  elec- 
trons are  not  considered.  This  simplifying  assumption  has  been  checked  in 
cylindrical  geometry  ns  in-  *ho  nVNACYl.  code  on  an  evacuat'd  cylinder  with 
differences  in  the  elec  . c i u-iJ  near  the  cylinder  axis  of  around  SOf . 

There  should  be  even  less  t'fect  in  the  spherical  case  because  the  shape 
of  the  emitting  surface  causes  less  charge  to  be  directed  toward  the  axis 
of  the  problem. 

To  have  the  capability  of  varying  the  zone  spacing  in  the  r coordinate, 
the  generalized  coordinate  q^  is  related  to  r hv 

r = f(qj)  , 

and  the  other  coordinates  are 

0 = q?  » 

* * is  • 
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The  restriction  to  two  dimensions  gives  the  line  element  as 
Qj  Aq^  * 2irr  sin  0 . 


In  the  other  direction,  the  line  element  is 


Ql  Aqj 


In  the  0 direction,  the  line  element  is 

q2  *q 2 = ffqj)  • 

A0 

For  the  particular  case  where  0-0,  we  take  3 a — . 

The  particle  motion  is  calculated  in  this  coordinate  scheme  using  tho 
relativistic  generalization  of  Newton's  law: 


where  P is  the  particle  momentum. 
P = y mv  , 


In  spherical  coordinates  with  azimuthal  symmetry,  the  force  equations 
become 

Yq 

qE  = mY  - — - m 
^ r r r 


and 


• mYr  Yy 
<*E0  = mY0  ♦ -7- 


where 

yr»  Yq  = yvr,  yvQ,  respectively  , 
q = particle  charge, 
m = particle  mass 
r = particle  radial  position, 
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vr’  ve  * Particle  velocities  in  r and  © directions,  respectively. 

Making  the  substitution  ? = yt  and  integrating  over  one  time  step,  assuming 
fields  arc  constant  over  the  step,  results  in 


and 


p Yfi  &t 

^r  1 ®0 

Y = Y + — - At  + 4 — 

r rQ  m 2 r0 


C Y Y At 

v - v , qE®  At  1 ^ 90 

e - ¥e0  m • 2 r0 


where  the  subscript  ”0"  indicates  the  value  of  the  quantity  at  the  earlier 

time,  v and  are  obtained  from 
r 0 

\ * V* 

and 

VG  = VY  * 


and  are  used  to  calculate  the  new  particle  position 


where 

r = (rQ  + r)/2  . 


This  completes  the  mathematical  description  of  the  DYNASPHERE  code. 


APPENDIX  C 


UPPER  LIMIT  FOR  SGEKP- INDUCED  SURFACE  CURRENT  ON  A SPHERE, 

IN  TERMS  OF  THE  EMISSION  CURRENT,  FOR  LONG- PULSE -WIDTH  EXCITATION 


The  surface  current  on  a sphere  is  related  to  the  sheet  current  X 
amp/m,  which  flows  on  the  sphere  surface.  Application  of  the  continuity 
equation  to  an  element  of  surface  area  yields 


3o 

3t 


7 


t 


a) 


where  a is  the  surface  charge  per  unit  area,  K is  the  sheet  current  (in 

2 

amp/m),  and  Jem  is  the  emission  current  density  (in  amp/m  ). 

Suppose  that  the  emission  current  is  varying  slowly,  such  that  sig- 
nificant variations  in  its  amplitude  occur  slowly  in  relation  to  the  time 
required  for  light  to  trivel  around  the  sphere.  In  this  case,  a will  be 

nearly  uniform  and  equal  to  the  net  charge  on  the  sphere  divided  by  the 

2 

sphere  surface  area  (4itR  ).  The  sheet  current  K will  be  a function  of 
position  on  the  sphere,  but  this  spatial  dependence  will  be  nearly  invar- 
iant in  time  for  a slowly  changing  emission  current.  For  these  conditions, 
Eq.  1 can  be  simply  integrated. 

We  define  the  surface  current  K about  a perimeter  2ttR  sin  0,  as  shown 
in  Figure  C-l.  Since  K(Q)  is  constant  about  this  perimeter. 


s(0)  * 2wR  sin  0 K(Q)  . (2) 

Integrating  Eq.  1 over  the  sphere  surface  between  0=0  and  0 = 0Q, 
one  easily  finds 

2wRZ  (1  - cos  eQ)  — 2^.  + K(0)  2itR  sin  0 = Iem(©0)  . (*) 

4iw 

where  Q is  equal  to  the  total  emission  current  and  I (9n)  is  the  integral 

cm  u ^ 

of  the  emission  current  density  between  0=0  and  6 * 0g.  Therefore,  using 
Eq.  2,  Eq.  3 becomes 

MV  * Weo>  [?  <>  * c°s  V]  • ("J 

Note  that  Ig  approaches  zero  for  9g  = 0 and  for  6g  = 180°,  as  required  by 
symmetry. 


97 


The  relations  of  Eq.  4 may  be  slowly  varying  (overall)  in  time,  but 
for  the  quasi-static  conditions  assumed,  Eq.  4 holds  during  any  short 
time  interval. 


Any  spatial  function  may  be  postulated  for  the  emission  function 
Iem(0).  If  Qq  is  large  enough  to  include  the  entire  photo-emitting  area, 


VV  " I 'em  " • >--os  V 


(5) 


For  half-sphere  emission,  Eq.  5 holds  for  8Q  > 90°. 

For  uniform  emission  (I  (9)  ■ const],  the  maximum  surface  current 

cm 

occurs  at  0 * 90°.  For  a cosine  spatial  distribution  of  the  emission  over 
a half-sphere,  the  maximum  surface  current  occurs  at  a slightly  smaller 
angle.  For  emission  over  less  than  a half-sphere  surface,  the  maximum  sur- 
face occurs  near  the  edge  of  the  emission  area,  and  can  approach  the  value 
of  the  emission  current  for  the  case  of  emission  from  a small  spot  near 
0 = 0. 

Returning  to  the  case  of  half-sphere  emission,  with  a uniform  spatial 
distribution  of  the  emission,  it  is  simple  to  show,  from  Eq.  4,  that 

VV  5 7 'em 

and  that  the  peak  value  occurs  at  Q = 90°. 


RT-09936 


Figure  C-l. 


Schematic  illustration  for  definition  of  surface 
current  in  amperes 
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APPENDIX  D 


DESCRIPTION  CF  THE  PHYSICS  AND  MODELING  USED  IN  THE  CALCULATION 
OF  THE  ELECTRIC  AND  MAGNETIC  FIELDS  IN  THE  DYNACYL  COMPUTER  CODE 


Passage  of  photon-induced  electrons  through  a gas  generates 
electromagnetic  fields,  ionizes  the  gas,  and  (through  the  action  of  the 
fields)  causes  currents  to  flow  in  this  partially  ionized  gas.  The  method 
of  calculating  these  fields  and  currents  is  described  here,  assuming  that 
the  motion  of  the  phot oo-^ene rated  electrons  (referred  to  as  primary 
electrons),  as  well  as  the  ionization  rate  Hue  to  these  eiectrons,  is 
known.  Of  course,  the  treatment  is  equally  valid  if  electrons  are  injected 
directly  into  the  region  of  interest,  in  which  case  these  injected  electrons 
are  called  primary  electrons. 

The  two  aspects  of  the  problem  considered  here  are  the  generation 
of  the  fields,  and  the  generation  and  behavior  of  the  ionized  gas.  The 
fields  are  calculated  for  a cylindrically  symmetric  region  enclosed  in  a 
perfectly  conducting  can.  The  laws  governing  the  generation  of  the  fields 
are  Maxwell's  equations.  Motion  of  the  primary  electrons  ionizes  the  back- 
ground gas  by  generating  electron  ion-pairs.  The  latter  electrons  are 
referred  to  as  secondary  electrons.  The  motion  of  the  secondary  electrons 
can  cause  further  ionization.  The  motion  of  the  secondary  electrons  is 
described  by  an  empirical  d’-ift  velocity  and  the  rate  of  ionization  by  these 
electrons  is  also  empirical. 

Maxwell's  equations  are  solved  in  cylindrical  coordinates. 

We  assume  rotational  symmetry  about  the  cylinder  axis.  The  rotational 
symmetry  is  tantamount  to  saying  that  all  derivatives  with  respect 
to  azimuthal  ,*ngle  are  zero.  However,  this  still  leaves  the  possibility 
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of  currents  in  the  azimuthal  direction  (toroidal  currents).  The  assump- 
tion that  there  is  no  preferred  direction  of  rotation  eliminates  this 
possibility.  The  fore  of  the  reduced  equations  is 
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These  are  the  equations  to  be  solved  for  the  radial  and  axial 
components  of  the  electric  field  (Er  and  E£)  and  the  azimuthal  magnetic 
field  (H^)  written  in  MKS  units.  The  medium  in  which  the  fields  are 
generated  is  such  that  the  medium  has  the  same  permittivity  and  permea- 
bility as  a vacuum.  The  region  of  interest  is  surrounded  by  a perfectly 
conducting  can.  It  is  assumed  that  all  fields  are  initially  zero,  so 
that  the  three  equations  above  imply  the  divergence  equations 


7 • €qE  * p 


7 • u H = 0 


if  the  charge  density,  p,  is  determined  by 


1 


From  the  form  of  the  equations,  it  is  clear  that  they  represent  the  time 
evolution  of  the  fields  driven  by  the  current.  Thus,  it  is  the  current  which 
must  be  specified  to  determine  the  fields.  The  current  is  due  to  the  motion 
•>f  the  electrons;  ions  are  considered  to  be  stationary  on  tinr  scales  of 
interest  here.  Those  electrons  which  are  moving  inertip'.ly  are  classified 


as  primary  electron?  and  those  whose  motion  is  better  described  as  drifting 
through  the  background  gas  are  classified  as  secondary  electrons. 

The  partially  ionized  background  gas  (air)  is  treated  as  a con- 
ductor. The  conductivity  is  equal  to  the  current  density  divided  by  the 
electric  field*  here  considered  e scalar.  Since  the  current  density  is 
the  charge  on  an  electron  multiplied  by  the  product  of  the  electron  dens  ty 
and  the  electron  drift  velocity,  the  latter  two  quantities  must  he  known. 

The  secondary  electron  density  must  be  calculated  by  integrating 
in  time  the  continuity  equation 


3N 
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where  No  is  the  secondary  electron  density,  v is  the  electron  velocity,  and 
S is  a source  term  (there  is  no  sink  for  the  electrons  considered  here). 

The  source  term  has  two  contributions.  The  first  is  the  ioniza- 
tion rate  due  to  the  primary  electrons  which  is  calculated  while  integrating 
the  equations  of  motion  of  the  primary  electrons.  The  second  contribution 
is  due  to  the  ionization  by  the  secondaries  themselves.  For  the  rate  coef- 
ficient for  this  last  process,  we  approximate  the  experimental  results  of 
Reference  1 by 


Rate  for  ionization  by  secondaries 


1.62  x 10* 13  m3/sec 


where  E is  the  magnitude  of  the  electric  field  and  N is  the  neutral  number 
density  in  MKS  units. 

Although  the  results  are  for  the  hydrogen  molecule,  in  the  region 
of  interest  there  should  not  be  considerable  difference  between  hydrogen 
molecules  and  nitrogen  molecules,  the  prime  constituent  of  air.  For  low 
E/N  ratios,  the  approximation  is  somewhat  low  compared  to  the  experimental 
results  for  nitrogen  molecules  given  in  Reference  2. 
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The  experimental  data  for  the  drift  velocity  are  presented  in 
Reference  2 , and  approximated  by  the  expression 


Drift  velocity  « 


),i  aJO^JE/NI 


♦ 8 x 10S  (E/H)  2 exp 


■1.7S  x 10 


‘it 


(»/»«<) 


litis  is  the  last  bit  of  information  needed  to  describe  the  background  gas. 
The  conductivity,  a,  is  thus  determined. 

The  motic.  of  the  primary  electrons  is  treated  separately  from  the 
pressure  effects  and  Maxwell's  equations.  This  separate  treatment  pro- 
vides both  the  primary  current,  and  the  rate  at  which  the  primary 
electrons  are  ionizing  the  gas  to  provide  the  secondary  electrons.  The 
primary  electrons  are  considered  as  high  energy  (greater  cnan  10  rydbergs) 
and,  thus,  the  Bethe  form  of  the  energy  loss  of  a particle  can  be  used 
(References  3 and  4)  to  derive  the  ioni  at ion  cross  section  if  one  assumes 
3 rydbergs  are  lost  per  ionization.  T’ e cross  section  for  nitrogen  mole- 
cule is 


o ionization 


where  Z is  the  ' »ic  number,  t/.ken  to  be  ? for  air,  h/n.±  is  the  Compton 
wavelength  (divided  by  2ir),  a .s  the  fine  structure  constant,  and  3 is  the 
electron  velocity  Jivided  by  the  speed  of  light.  The  expression  is  non- 
relativistic.  Of  course,  the  above  cross  section  can  be  used  to  calculate 
the  rate  of  slowing  of  the  primary  electrons. 

The  total  current  is  then 


J * J ♦ oE 
P 

This  fully  determines  the  solution  and  it  remains  to  describe 
the  treatment  of  the  numerical  procedure  rnd  the  botaidary  conditions. 
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The  centering  in  tine  end  apace  it  summarized  below. 


Coordinate 


s r t 

b c b 
c c c 
b b c 
c c b 
b b b 
c c c 


The  tern  b indicates  evaluative  at  a boundary  and  c indicates  evaluation 
at  a center  in  the  appropriate  coordinate.  This  cratering  scheme  is  some- 
what unfortunate  in  that  it  lacks  symmetry  in  the  coordinates  and  causes 
undue  difficult/  in  the  boundary  conditions.  However,  these  difficulties 
have  been  oveicoaie.  Actually,  the  above  centering  is  not  precisely  true 
for  the  current  since  the  conduction  current  is  calculated  using  the  new 
value  of  the  electric  field.  This  implicit  scheme  insures  numerical 
stability. 

For  a perfectly  conducting  can,  the  boundary  conditions  require 
that  the  electric  fields  arc  normal  to  the  can  the  magnetic  field  is 
parallel  to  the  can  at  the  boundaries.  The  electric  field  in  the  z direc- 
tion is  thus  zero  at  the  radius  of  the  can.  However,  due  to  the  centering 
used  here,  the  radial  electric  field  is  not  evaluated  on  the  top  and  bot- 
tom of  the  can,  but  is  evaluated  half  a zone  inside  the  can.  Thus,  a 
field  equal  and  opposite  is  assumed  to  exist  half  a zone  outside  the  can, 
and  this  is  used  to  calculate  tho  time  derivative  of  the  magnetic  field 
at  the  ends  of  the  can.  The  fact  that  H is  the  only  nonzero  magnetic 
field  automatically  fulfills  the  boundary  condition  on  the  magnetic  field. 

A further  difficulty  occurs  at  the  axis  where  the  radius  is  zero. 
To  evaluate  the  curl  of  the  magnetic  field  there,  resort  must  be  made  to 
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the  definition  f the  curl  which  is  the  line  integral  a round  the  boundary 
of  the  surface  divided  by  the  area  of  the  surface . this  provides  the 
necessary  quantity. 

Surface  currents  and  charges  «tre  secondary  quantities  in  that 
they  are  derived  fro*  other  fields  and  currents.  Surface  currents  are 
calculated  by  requiring  then  to  be  such  as  to  make  the  magnetic  fields 
tero  outside  the  can. 

This  completes  the  seminary  of  the  treatment  of  the  background 
gas  for  pressure  effects  and  the  calculation  of  the  elect rtmagnetic  fields, 
fhe  numerical  treatment  is  straightforward  and  stable  as  long  as  the  time 
step  is  chosen  so  that  light  csnnoi  travel  across  more  than  half  a a one  in 
one  time  step. 
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DESCRIPTION  OP  PHYSICS  AND  M)Wl.tN3  USED  IN  I1EIBS 
AND  CURRENTS  IN  HIE  TSFHERIi  CONTUTER  CODE 


This  appendix  describes  the  physics  and  Modeling  employed  in  the 
TSHtEKE  computer  code  to  obtain  charge  distrihut  ions,  electric  fields,  and 
surface  currents.  Charge  distributions  are  obtained  by  following  individual 
"particles**  of  charge  which  were  emitted  by  the  inner  sphere.  Electric 
fields  are  calculated  from  the  charge  distribution  employing  a Green’s 
function  which  is  integrated  over  the  distribution  each  time  step.  Surface 
currents  are  obtained  from  the  electric  field  at  the  inner  sphere  boundary 
employing  th-.s  continuity  equation  which  relates  the  spatial  gradient  of  the 
surface  current  to  the  time  rate  of  change  of  the  normal  electric  field  at 
the  surface.  Each  of  these  methods  is  described  in  detail  in  the  following. 

The  charge  density  is  obtained  in  TSPHERE  by  following  individual 
particles  of  charge.  These  particles  are  acted  on  by  the  electric  fields 
through  the  appropriate  force  equations.  Their  positions  are  updated  each 
time  step.  The  particles  are  then  collected  into  the  different  spatial 
cells  to  obtain  the  spatial  distribution  of  the  charge  density.  The  cal- 
culation is  self-consistent  in  that  the  particles  move  consistently  in 
fields  which  were  calculated  from  the  charge  distribution  of  the  previous 
time  step. 

Relativistic  effects  are  included  in  the  particle  position  updating 
scheme  used  in  TSPHERE,  The  equations  are  obtained  by  integrating  the 


equations 


where  P is  the  particle  momentum. 


ftniig  pm  Mat 


- , .1  ^ 


'•1,  - ;v  '.-  •£ 


&i  "■"  ■ ■ ■ ■ - • 


P » y n«v 


A - *)' 


In  spherical  coordina  es  with  azimuthal  symmetry,  the  force  equations  become 


qEr  = mvr  - ~ m 


, mv  va 

qE0  = mv6  + r"~ 


where 


q = particle  charge  , 

m = particle  mass, 

r = particle  radial  position, 


vr»  vg  s particle  velocities  in  r and  0 directions,  respectively. 


» ®"r  y 

Making  the  substitution  Y = yv  and  integrating  over  one  time  step,  assuming 
fields  are  constant  over  the  step,  results  in 


qEr  1 Y0O  At 

Yr  = Yr~  + — - A,  + ~ — 

0 m 2 r0 


P Y Yq  At 

v . , . qEe  1 ro  eo 

0 0Q  m 2 rQ 


where  the  subscript  "0"  indicates  the  value  of  the  quantity  at  the  earlier 


time.  v and  vQ  are  obtained  from 
r y 


vr  = Vr.Y 


v = Y /y 

e e/Y  • 


and  are  used  to  calculate  the  new  particle  position 

r - r0  * f(Vro  * Vr)M  • 

0 ■ 0o  * 7 (V0o  * V0)  at/7  • 

where 

r - (rQ  + r)/2 

Electric  fields  are  obtained  in  TSPHERE  by  solving  Poisson's 
equation  each  time  step: 

V2  $ = - p/t  (mks  units)  (E.l) 

where  $ - electric  potential 
p = charge  density 

The  Green's  function  technique  is  used  for  the  solution.  The  space  between 
the  spheres  is  broken  into  radial  and  angular  zones.  Each  resulting  cell 
is  actually  a donut-like  shape  when  it  is  rotated  about  the  axis  of  symmetry. 

is  obtained  analytically  from  Eq.  E.l  for  individual  rings  of  infinitesimal 
thickness  at  the  locations  of  every  zone  in  the  mesh  by  solving  the  equation 

V-  * = - 6(r  - c)  6(9  - a)  . (E.2) 

The  terms  r,9  is  the  field  point  and  c,a  is  the  source  point. 

The  solution  to  Eq.  E.2  is  in  terms  of  an  infinite  series.  The 
series  is  differentiated  in  each  direction  term  by  term  to  obtain  expressions 
for  the  electric  fields.  The  resulting  expressions  are: 
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In  the  above  expressions,  the  following  definitions  apply. 


a 


I, 

u 

l- 

E-  ' 
r 

iV 


GrijU,  G0i j kfi.  = 


e. 

j 


p 

n 


Green's  functions  relating  source  points 
c^,  cx^,  to  field  points  r^ , 0^  , for  radial 
and  electric  fields,  respectively  angular 
field  points 
source  points 

inner  and  outer  sphere  radii 
Legendre  polynomial  of  order  n 


These  expressions  are  then  integrated  term  by  term  over  each  zone  volume 
to  obtain  Green's  functions  consistent  with  our  particle-in-cell  ,*-thod  for 
obtaining  the  charge  distribution.  The  electric  fields  are  then  obtained 
from  the  Green's  functions  by 


no 


E-44  * ]C  Grijki  p 


rij  Z*  rijki  ‘'ici  . 2 

k,i  **e0  rl 


where 


E0ij  " 2 G0ijkA  Pkt 


Erij  * E0ij 


radial  and  axial  electric  fields  at  field 
point  r^, 

charge  density  at  source  point  c. , a 


k*  a 


net  charge  on  inner  sphere 


The  surface  current  on  the  inner  sphere  is  obtained  from  the  net 
current  and  normal  electric  field  there.  It  is  obtained  by  integrating 
Ampere's  law: 


V x H = 


t ^ 3E 

J * Eo  n 


both  in  space  and  time.  Remembering  that  the  tangential  component  of  H at 
the  boundary  is;  equal  to  the  skin  current  there  (in  units  of  current  per 
unit  length),  we  obtain: 


I = 2ttR2 


JsinOd0 


sin0dO 


Refer  to  Figure  3-3  in  Section  3 for  definitions  of  the  coordinates. 

J is  the  net  current  at  the  position  0 on  the  inner  sphere  Er  is  the  normal 
electric  field  there.  I is  the  surface  current  at  0 in  amps. 

In  the  code,  the  quantities  for  the  net  current  J and  the  time 
derivative  of  the  normal  electric  field  are  averaged  over  a number 
of  time  steps  in  order  to  produce  smoothly  varying  surface  currents.  This 
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step  is  necessary  at  higher  fluences  because  of  the  particle  nature  of  the 
code.  In  a given  time  step,  many  particles  might  pass  through  a small  area 
about  the  position,  0,  The  next  step  there  might  be  very  few.  The  result 
is  thnc  ’.he  quantity  J can  oscillate  in  time.  At  high  fluences,  the  two 
terms  in  "he  integral  for  I are  large  and  their  sum  small  so  large  oscillations 
in  either  term  can  render  the  sum  completely  unless  appropriate  steps  such 
as  just  described  are  taken. 


